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16094 LOCAL DESTRATIFICATION LAKES 


KEY WORDS: Dispersion; Hydraulics; 
Hydrodynamics; Lakes; Mixing; Models; Reservoirs; Stratification; Water 
resources 


ABSTRACT: numerical simulation the flow field the jet induced axial- 
flow propeller pump performed solving the governing equations the flow field. 
The two-dimensional Los Alamos SOLA prediction technique has been expanded 
include diffusion well buoyancy forces axisymmetric stratified flow field. 
Comparison with hydraulic models shows that this numerical simulation the local 
destratification phenomenon useful for the prediction jet penetration depth 
the dilution factor (released water quality). The main dynamic effects are modeled 
adequately show the same trends the available physical data. The prediction 
procedure represents low-cost basic tool show the influence design parameters 
the flow field. empirical expresssion for prediction the penetration depth 
suggested the basis experimental and computational results. 


REFERENCE: Busnaina, Ahmed A., Lilley, David G., and Moretti, Peter M., 
“Prediction Local Destratification Lakes,” Journal the Hydraulics Division, 
ASCE, Vol. 107, No. HY3, Proc. Paper 16094, March, 1981, pp. 259-272 


16099 EROSION RECTANGULAR LONG CONTRACTION 


KEY WORDS: Bed load; Channel design; Channel erosion; Channels 
(waterways); Contraction; Erosion; Friction; River engineering; Scour 


ABSTRACT: experimental investigation into the general erosion within long 
contraction problem reported. The experimental results are used verify Straub’s 
theoretical one-dimensional model for general erosion. Local scour near the side walls 
the inlet region was also measured and the results are reported. Two different sizes 
sand were used the experiments. 


REFERENCE: Gill, Mohammad Akram, “Bed Erosion Rectangular Long 
Contraction,” Journal the Hydraulics Division, ASCE, Vol. 107, No. HY3, Proc. 
Paper 16099, March, 1981, pp. 273-284 


16110 MODEL FOR WIND INDUCED FLOW 


KEY WORDS: Channel flow; Depth; Eddy currents; Models; Numerical 
analysis; Surface waters; Three dimensional flow; Wind; Wind forces 


ABSTRACT: Based compilation current data the surface, the classical 
Ekman theory inappropriate for calculating wind generated currents near the surface. 
The cause found the assumption Ekman vertical.eddy viscosity, Ny, 
which constant over depth. numerical model described briefly, which allows 
for variation over depth, well the and directions. The model was 
developed using Galerkin technique. The model can include sharp, wind driven 
profile near the surface. The method computationally reasonable, requiring about the 
same number computations vertically averaged formulation. relationship for 
choosing versus depth presented. Experimental data are compared the 
numerical model results. 


REFERENCE: Pearce, Bryan R., and Cooper, Cortis K., “Numerical Circulation 
Model for Wind Flow,” Journal the Hydraulics Division, ASCE, Vol. 107, 
No. HY3, Proc. Paper 16110, March, 1981, pp. 285-302 


16103 CELEBRATION KARAJI’S HYDROLOGY 


KEY WORDS: Ground water; Groundwater hydrology; History; Hydrologic 
cycle; Underground conduits; Underground structures 


ABSTRACT: brief account Mohammed Karaji, Persian scholar the tenth 
century, and his contributions science, especially ground-water hydrology, 
presented. Karaji wrote the first algebra text presented with symbolism. engineering, 
invented instruments that proved useful the surveying and the tunneling 
underground chanels for conveying water. Most importantly, expounded the basic 
principles hydrology, understanding the proportionality law between the force 
gravity and the mass body. Karaji noted that ground-waters, like surface waters, 
can stagnant, and that still waters can located under vast deserts and lowlands. 
recently discovered text Karaji “hidden the oldest known text 
ground-water hydrology. 


REFERENCE: Pazwash, Hormoz, and Mavrigian, Gus, Celebration 
Hydrology,” Journal the Hydraulics Division, ASCE, Vol. 107, No. HY3, 
Proc. Paper 16103, March, 1981, pp. 303-309 


16102 TURBULENCE MEASUREMENT STUDY 


KEY WORDS: Electromagnetism; Flow measurement; 
Hydromechanics; Instrument sensitivity; 
Turbulence; Turbulent flow; Water current; Water metering 


ABSTRACT: The capability modified commercially available electromagnetic 
water current meter was evaluated for measurement parameters necessary 
quantify turbulent flow. Values turbulence intensity were measured, and Fourier 
analyses their oscillograph recordings were performed obtain the spectral 
distribution the energy. Field tests were then conducted three locations selected 
the basis their differences flow regime. The results showed that values 
turbulent intensity for flows between very high light turbulence levels ranged from 
0.71 0.06, respectively, and the Fourier analysis the data shows the energy 
concentrated generally below Hz. This instrument can used obtain relative 
values turbulent intensity and spectral information from field measurement. 


REFERENCE: Griffith, Franklyn, III, and Grimwood, Charles, 
Measurement Study,” Journal the Hydraulics Division, ASCE, Vol. 107, No. HY3, 
Proc. Paper 16102, March, 1981, pp. 311-326 


16144 ULTIMATE DIMENSIONS LOCAL SCOUR 


KEY WORDS: Analysis (mathematics); Comparative studies; Erosion; 
Hyperbolic functions; Linear functions; Logarithms; Scour; Time 
measurement; Urban highways 


ABSTRACT: mathematical method presented for determining the ultimate 
dimensions local scour from relatively short term measurement the progression 
scour with time. The method involves plotting the logarithm the rate change 
the scour dimension against the logarithm time. hyperbolic curve fitted the 
plotted data, and the asymptote the hyperbola used determine the ultimate 
scour dimension. example used compare the hyperbolic logarithmic method 
with the linear logarithmic and linear semilogarithmic methods (in which the scour 
depth increases without limit) which require estimate the time required reach 
practical equilibrium. The hyperbolic logarithmic method predicts both the time 
progression scour and the ultimate scour dimensions. 


REFERENCE: Blaisdell, Fred W., Anderson, Clayton L., and Hebaus, George G., 
“Ultimate Dimensions Local Scour,” Journal the Hydraulics Division, ASCE, 
Vol. 107, No. HY3, Proc. Paper 16144, March, 1981, pp. 327-337 


16124 THREE-PARAMETER PROBABILITY DISTRIBUTIONS 


KEY WORDS: Flood frequency; Floods; Frequency distribution; Hydrology; 
Low flow; Probability density functions; Probability distribution functions; 
Probability theory; Statistical analysis; Weibull density functions 


ABSTRACT: The lognormal, Weibull, Pearson type and log Pearson type each 
three-parameter distribution, were evaluated generalized fashion terms the 
dimensionless variate (K=X/vx, which equals randon variable, and equals 
its mean) which has population mean unity. The bounds the distributions, areas 
the portions distributions the negative region variate when they enter such 
regions, and the differences some important quantiles among the four distributions, 
are presented. The four distributions become less applicable for hydrologic frequency 
analysis they deviate more and more from their two-parameter counterparts 
(lognormal the case log Pearson). When they have well-applicable properties, 
their quantile values differ little for some all distributions indicating that choice 
distribution makes little difference. Some guidelines are provided for selecting the best 
applicable distribution for given hydrologic sample. 


REFERENCE: Rao, Probability 
Distributions,” Journal the Hydraulics Division, ASCE, Vol. 107, No. HY3, Proc. 
Paper 16124, March, 1981, pp. 339-358 
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Thermal stratification may occur water reservoirs during the hot months 
summer, and then three main layers water with different characteristics 
may observed. The epilimnion, the top layer, contains warm low-density 
water (usually rich oxygen because atmosphere reaeration and photosynthe- 
sis, and thus considered high-quality water). The hypolimnion, the bottom layer, 
consists cold high-density water (often poor oxygen and thus considered 
low-quality water). The region rapid temperature change (between the other 
two layers) called the thermocline metalimnion. This stratification presents 
serious problem for reservoirs with low-level release structures that the 
quality water released, characterized its oxygen content, may poor, 
since most the exit flow comes from the bottom layer water the reservoir. 
Since many old reservoirs have release structures located near the bottom, 
there indeed problem finding suitable and economically feasible alleviation 
techniques. 

Artificial destratification, structural modification, and localized mixing are 
three possible alleviation techniques that can used improve the release 
water quality. Artificial destratification can either mechanical pumping (with 
assorted piping) diffused-air pumping. These mixing devices, however, require 
substantial amount energy destratify large body water, (see Ref. 
13). Structural modification the dam involves elevating the release gate position 
order allow some the water drawn from the epilimnion. This 
method, although effective, extremely costly. the other hand, localized 
mixing epilimnion water into the hypolimnion has been found effective 

Student, School Mech. and Aerospace Engrg., Oklahoma State Univ., Stillwater, 
Okla. 74074. 
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and economical enhancing the quality water release from low-level release 
gates. There has been continued interest local mechanical destratification 
reservoirs improve water quality. low-energy axial flow propeller may 
positioned just below the water surface provide downward-directed 
jet water and thereby locally mix the reservoir the vicinity the release 
structure the dam. The high-quality epilimnion water transported downwards, 
and the effectiveness pump this type has been illustrated field 
experiments Pine Creek Reservoir, Okla. (8), and Lake Okatibbe, Miss. (6,7). 
Fig. shows schematic the mixing produced when there exit beneath 
the propeller. The jet plume penetrates some distance below the level 
the thermocline. There mixing epilimnetic and hypolimnetic water the 
shear layers the plume, and the resulting intermediate-density water spreads 
out lens-like gravity flow the level the thermocline. 

Fig. shows the practical application localized mixing the proximity 
the release structure dam with low-level release gate. The flow field 
produced fully three-dimensional. However, axisymmetric two-dimensional 
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FIG. 1.—Schematic Typical Pump and Flow Field Produced without Exit 
Flow 


simulation this phenomenon shown Fig. this, the release flow 
represented circular opening the bottom the flow field, directly 
beneath the propeller. This does not represent correctly the position the 
exit, since practice off-axis position, but provides useful 
simplification the true problem. either case, the release flow can entrain 
some all the jet plume water from the surface. the arrangement 
successful, the release flow will consist mostly epilimnetic water, along 
with some hypolimnetic water that enters the plume through mixing, and the 
dilution factor (as defined later) will near one. the velocity the downward 
jet insufficient for penetrating deeply into the stratification, only hypolimnetic 
water will withdrawn, and the dilution factor will zero. the jet velocity 
adequate, but the propeller flow rate small compared the release flow 
rate, the plume located too far from the release gate, then the dilution factor 
will have some intermediate value. 

Hydraulic model experiments have enabled considerable success achieved 
understanding the phenomena (11). Recent numerical simulations and 


LOCAL DESTRATIFICATION 


FIG. Typical Propeller Pump and Flow Field Produced with Exit 


Flow via Low-Level Release Gate: (a) Release Structure Dam; (b) Directly Beneath 


local mixing supplement the experimental program, and demonstrate substantial 
promise for improvement modeling the transient mixing behavior during 
destratification. Advantages being able calculate (rather than experiment 
the laboratory the field) are enormous terms their cost/benefit 
ratio. Such procedures confirm and extend the available experimental data 
larger set cases, and supplement the rather limited theoretical treatment 
associated problems (see, e.g., and 5). 


Several fundamental questions arise the correct modeling the physical 
phenomena. First, consider the question modeling turbulence. Over the range 
parameters studied, flow visualizations the laboratory suggest that the 
controlling factor for the dilution the release water will epilimnetic water 
the entrainment the downward plume into the release flow, rather than 
the mixing process the vicinity the plume (which the governing mechanism 
for whole-lake mixing). Thus, the modeling jet penetration depth the 
one hand, and critical release flow rates the other, are key items the 
representation the physical phenomena. The viscous turbulent mixing 
processes the shear layers are less important, and not greatly affect the 
jet velocity which the dilution factor increases from low high value. 
the numerical model, this has been confirmed varying the viscosity from 
laminar value constant large viscosity with minor effect the results. 
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That such value appropriate for round turbulent free jet well-known 
(12). Since the computed effects were minor, may concluded that precise 
modeling the turbulence not crucial obtaining good results. The present 
work incorporates laminar value for the viscosity, coupled with turbulent-flow 
value for the Schmidt number (equal unity) used the species equation. 
The result good similitude between temperature-stratified field tests, chemi- 
cally-stratified laboratory tests, and numerical predictions. 

Another question the difference between axisymmetric problems which 
the present theory based, and asymmetric problems which the experimental 
data Moon (10) were obtained, and with which results are compared. Identical 
results cannot expected, but the work reported shows that similar parametric 
trends are observed both geometries. This would undoubtly not the case 
the geometry the asymmetric case were longer compact, with the release 
gate laterally far away from the plume. 

The numerical simulation the flow field the jet, induced axial-flow 
propeller pump, performed solving the governing equations the flow 
field, i.e., the equations conservation mass (continuity), momentum (in 
axial and radial directions), and species diffusion. Density changes are slight 
and allow the flow field simulated incompressible flow with the 
addition upward buoyancy force. The equations cylindrical (axisymmetric) 
coordinates may taken conservative form as: 


2 


am, m,U 


which the coordinates the radial and axial directions, respectively; 
pressure deviation from hydrostatic pressure; dynamic viscosity; 
gravitational acceleration; mass fractions fluids and respectively; 


— + 
am, 
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the Schmidt number. 

With these equations, vertically downward-directed jet the stratified 
environment can simulated, and proper account taken species mixing 
and buoyancy. The computational method based the transient two-dimen- 
sional Los Alamos SOLA prediction technique, described Ref. which 
finite difference scheme based the Marker and Cell method. The finite 
difference equations are based directly the primitive pressure-velocity vari- 
ables, rather than the associated stream function and vorticity variables, 
and this proves particularly convenient the imposition boundary 
conditions. 

The flow domain shown Fig. represents the physical problem. has 
vertical axis symmetry and contains downward-flowing jet fluid from 
the propeller. Initially, two fluids occupy positions above and below the interface, 


FREE SURFACE 


PROPELLER 


INTERFA 
AXIS 
SYMMETRY 


x 


FIG. 3.—Axisymmetric Flow Domain 


the height the interface), and vice versa. Fig. shows how the cylindrical 
region divided into equal-sized rectangular cell divisions. This solution domain 
complemented layer cells all sides, allow easy simulation 
the required boundary conditions. single cell this mesh enlarged 
Fig. which shows the location each field variables, and other 
variables this cell. The pressure, and mass fraction fluid m,, are 
located the center each cell and the radial and axial velocities lie directly 
the physical boundaries the solution domain, while the pressure and mass 
fractions are displaced half cell interval inside the flow field. this way, 
the exterior fictitious are particularly convenient when applying the boundary 
conditions. 

Fig. also shows how inlet and exit flows are handled. The exit release 
flow via circular opening directly beneath the propeller, with the flow 
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rate specified priori. The available volume the flow domain very limited. 
allow the outlet release water without dramatically decreasing the fluid 
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FIG. Finite Difference Variables Typical Cell 


level, there compensating incoming flow which distributed among the 
top and bottom layers amounts equal the epilimnetic and hypolimnetic 


PROPELLER 

FIG. 4.—General Mesh Arrangement 
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water released. The peripheral inflow two different elevations such 
way that epilimnetic water enters through circumferential opening into the 
upper layer and does not disturb the density profile. The large circumferential 
opening allows the inflow low velocity order avoid any disturbance 
the flow field. Similarly, hypolimnetic water fed into the bottom layer. 

The solution algorithm consists time-march procedure. each step 
the march, velocity guesses are calculated explicitly from their respective 
conservation momentum equations. Application the diffusion equation 
(conservation chemical species) the low-density fluid (initially the top 
layer fluid the solution domain), enables the respective proportions 
fluid and fluid (the high-density fluid initially forming the bottom layer 
fluid the solution domain), calculated, and thus, the density 
all locations. The newly-calculated velocities are considered first approx- 
imations the new time level, since the continuity requirement has not yet 
been imposed. The approximations are improved iterative fashion, adjusting 
cell pressures and velocities satisfy the continuity equation. When convergence 
achieved, current values pressure and velocity are accepted and taken 
initial values for the next time step. Associated useful calculations may now 
deduced terms these values. 

Boundary conditions solid boundaries are imposed the velocities and 
species concentration after each time step and after each pass through the mesh 
during pressure iterations. Convergence the steady-state solution established 
taking many forward time steps. The choice time increment must 
restricted (for stability) two ways. First, fluid should not pass through more 
than one cell one time step. Therefore must less than (usually 0.25 
times) the minimum cell transit time taken over all cells. Secondly, 
when nonzero value kinematic viscosity used, the momentum should 
not diffuse more than one cell one time step. Accuracy enhanced 
using small space and time intervals, the choice being determined 
hoc manner suitable trade-off between accuracy and economy computation. 
full examination the general finite difference technique and its application 
this problem are readily available (3,9). 


Comparison results with hydraulic laboratory model data Moon (10) 


shows that the numerical simulation local destratification phenomena 
useful predictive tool. The dilution factor (release water quality) follows: 


which the density; and subscripts and refer release, hypolimnion 
(bottom), and epilimnion (top) locations, respectively. The dilution factor 
found strong function the following. 


The densimetric Froude number: 


V. 
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which the initial jet velocity from the propeller; the gravitational 
acceleration; and Ap/p, the stratification. 
The nondimensional metalimnion (interface) location: 


which the metalimnion location measured from the bottom. 
The nondimensional flow rate: 


which represents the flow rate; and subscripts and refer the propeller 
and release locations, respectively. 
The nondimensional propeller diameter: 


which the propeller diameter. 
The nondimensional propeller depth: 


which propeller depth below the surface. 

standard base case taken with the following parameters (typical 
0.211, from which certain changes are made and individually assessed the 
subsequent analysis. Figs. 6-10 show the computed effect variation each 
these parameters the dilution factor, DF, while the other parameters 
maintained their respective base values. Appropriate experimental results 
are also shown where and observed that the main dynamic 
effects are modeled adequately and show the same trends the available physical 
data. 

Fig. the dilution factor, DF, shown strong function the 
densimetric Froude number, both predictive and laboratory data. The 
general result that release water quality (the fraction epilimnetic top 
water the exit stream) increases with higher values the Froude number, 
such conditions being achieved with higher jet velocities from the propeller, 
lower degree stratification the lake, both. The favorable comparison 
experiment and theory clearly seen. From the designer’s view, the important 
result the Froude number which the dilution factor rises sharply, since 
that determines the velocity required the propeller for the given conditions. 
This Froude number predicted very well the computer program. second 
point interest the maximal dilution factor achievable, which predicted 
within the scatter the experiments. 

The relation between the dilution factor, DF, and the normalized metalimnion 
location, presented Fig. The computation carried out for the 
base case which the relatively flat part the dilution factor 


(9) 
Q, 


HY3 LOCAL DESTRATIFICATION 267 


versus the Froude number curve (see Fig. 6), and results dilution factors 
having nearly maximal values. The location the metalimnion seen Fig. 
have strong effect the penetration depth the jet water from 
the propeller, and thus has effect also the ultimate dilution factor. 
the metalimnion (interface) height (above the base) increases, buoyancy forces 
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FIG. 6.—Dilution Factor, DF, Function Densimetric Froude Number, 
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oO 


FIG. 7.—Dilution Factor, DF, Function Metalimnion Location 


increase, and this will hinder the jet penetration, increase the dispersion 
the plume, and thus reduce the dilution factor. experimental data are currently 
available directly substantiate this predicted result. 

Fig. shows that the dilution factor, DF, strong function the flow-rate 
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Q*. Generally, release water quality improves with increase the 
value Q*, condition achieved with low release rates and high propeller 
flow rates. The figure shows the relationship found both the hydraulic laboratory 
model and computer simulation, and the agreement extremely good. Contrary 
first intuition, the dilution factor can exceed the flow-rate ratio. This illustrates 
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FIG. 8.—Dilution Factor, DF, Function Nondimensional Flow Rate, 
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FIG. 9.—Dilution Factor, DF, Function Nondimensional Propeller Diameter, 


the complex interaction between the propeller flow field and the release flow 
field, and shows the need for the application numerical computations. 

The effect changing the propeller diameter, but retaining the propeller 
flow rate, Q*, the same, shown Fig. Again, there good agreement 
between predicted and experimental evidence the effect normalized diameter, 
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D*, dilution factor, DF. These and other results show that values 
less than 0.1 are ineffective, and that value about 0.13 gives the best 
results. view the differences geometry, the discrepancies between 
hydraulic and numerical models are surprisingly small. 

The final parameter influence considered that depth the propeller 
below the surface. Fig. shows predicted, almost linear slight decrease 
dilution factor, DF, the nondimensional propeller depth, L*, increases, 
illustrating that the propeller lowered toward the metalimnion, less epilimnetic 
water driven the exit. Predicted results with values less than 
are unrealistic, since the numerical model constant surface location 
specified. practice, surface effects would occur when the propeller too 
high the water, with entrained vortex air bubbles being pumped downward, 
thus decreasing the water flow rate, jet penetration, and dilution factor. Experi- 
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FIG. 10.—Dilution Factor, DF, Function Nondimensional Propeller Depth 


mental data are not available for full range depths, but the figure includes 
some experimental data obtained with 1.0; 0.75; and 0.183. 
Though these are different parameters from the ones used the prediction, 
effects the proximity free surface are evident, and the predictions 
these cases are clearly not representative the physical situation. 

When there exit flow, (Q* «), the nondimensionalized jet penetration 


which the jet penetration depth measured from the surface; and 
the total depth, special interest. For example, then one 
assured that epilimnion water transported the lowest level, and release 
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gate placed there will enhance the epilimnion water transportation. interesting 
investigate the value function densimetric Froude number, 

empirical expression for prediction the penetration depth suggested 
the basis experimental and computation results (3), and deduced 

2 h 

which the jet downward velocity; the local jet center line downward 
velocity distance, below the propeller diameter, the local 
density; and empirical constant evaluated 3.5 from the experimental 
data Moon (10). That this equation exhibits good qualitative predictions when 
compared with available data seen Fig. 11, which shows normalized 
penetration depth (whose corresponding depth, below the propeller, may 
obtained setting equal zero the aforementioned expression) 
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FIG. 11.—Nondimensional Penetration Depth, Function Densimetric Froude 


function densimetric Froude number, for particular values dimensional 
propeller diameter, D*, initial metalimnion (interface) location, and norma- 
lized propeller depth, L*. expected, increases, also does the 
nondimensional penetration depth, Beyond Froude number unity, 
penetration the jet the floor the lake predicted. Similar results are 
also obtained with different values the parameters. 


The prediction procedure represents low-cost basic tool show the influence 
design parameters the flow field. Comparison with hydraulic models shows 
that this numerical simulation the local destratification phenomenon useful 
for the prediction the dilution factor (released water quality) and jet penetration 
depth. The main dynamic effects are modeled adequately show the same 
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trends the available physical data. The dilution factor and the nondimensional 
depth are functions the densimetric Froude number, metalimnion location, 
propeller flow rate relative the release flow rate, propeller diameter, and 
propeller depth below the surface. empirical expression for prediction 
the penetration depth suggested the basis experimental and computational 
results, and correlates well the available data within its range applicability. 
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The following symbols are used this paper: 
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constant, with the deduced value 3.5; 

propeller diameter; 

D,/H nondimensional propeller diameter; 

densimetric Froude number; 
gravitational acceleration; 

total depth model; 

depth below propeller; 

propeller depth below the surface; 

L/H nondimensional propeller depth; 

mass fraction fluid density, (epilimnion); 
mass fraction fluid density, (hypolimnion); 
pressure deviation from hydrostatic pressure; 
propeller flow rate; 

release flow rate; 

nondimensional flow rate; 

velocity radial direction; 

velocity axial direction; 

initial jet velocity; 

radial direction; 

axial direction; 

penetration depth (measured from the surface); 
Z,/H nondimensional penetration depth; 
metalimnion location (measured from the bottom); 
nondimensional metalimnion location; 
dynamic viscosity; 

initial density profile function 
epilimnion top water density; 

hypolimnion bottom water density; and 
Schmidt number. 
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BED EROSION RECTANGULAR LONG CONTRACTION 
Mohammad Akram ASCE 


When the width rectangular channel reduced building long side 
walls that the width within the side walls remains constant, the contracted 
channel may called rectangular long contraction. river channel within 
the parallel guide banks, built diversion works such weirs and barrages, 
practical example rectangular long contraction. The width river channels 
greatly reduced within the guide banks order economize the cost 
constructing diversion works. For example, the site the Mailsi Siphon, 
Pakistan (see Fig. 1), the natural width the Sutlej river reduced from 

The effect reducing the channel width increase the bed shear stress 
the flow, which responsible for causing considerable erosion within the 
contracted reach. order guard against the destructive effects excessive 
erosion within the reach, important able estimate the depth 
maximum erosion. The results laboratory study are presented here, which 
may hopefully used for the purpose practical design. 

One the earliest works which presented simplified one-dimensional theory 
the equilibrium depth erosion for rectangular long contraction due 
Straub (5). His work was later followed Ashida (1) and Komura (3). 
the Indo-Pakistani river engineering practice, suitable scaling factors (>1.0) 
are used together with the regime flow depth for estimating the depth maximum 
erosion long contractions. Essentials Straub’s one-dimensional theory are 
presented here and the final result will used fit the experimental data. 


Dept. Civ. Engrg., Ahmadu Bello Univ., Zaria, Nigeria. 

Note.—Discussion open until August 1981. extend the closing date one month, 
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One-Dimensionat Erosion 


With reference Fig. the equation continuity for water flow written 


which total discharge; channel width; average cross-sectional 
velocity; average flow depth; and superscripts and refer the upstream 
uncontracted and contracted channel reaches, respectively. assumed that 
the hydraulic resistance, both the upstream uncontracted and contracted 


with 
f pitch stone on siopes 
‘ @nd rubble apron 
toe 


Permanent earthworks 


FIG. View Mailsi Siphon, Pakistan; Note Layout Guide Banks Which 
Similar That Long Contraction 


FIG. Diagram Flow Through Long Contraction 


Prot of Temporary spur 
— 
~ 
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RECTANGULAR LONG CONTRACTION 
reaches, described the Manning formula, written follows 


which Manning’s coefficient bed friction; and energy gradient. 
The equation continuity sediment discharge written 


which unit sediment discharge. The rate sediment transport assumed 
given 


acceleration; specific gravity the bed material; average size 
bed material; bed shear stress; and critical shear stress the bed 
material. Assuming that any depression water level within the contraction 


negligible, and also ignoring the effect sorting, Eqs. 1-4 can solved 
give 


Eq. holds for flows which the sand bed the upstream uncontracted 
channel moving. For flows which the upstream uncontracted sand bed 


d, B, 6/7 7, 3/7 


The maximum depth erosion predicted the threshold movement 
the upstream bed, i.e., when 1.0, and given 


The highest value the known empirical formulas the 
Einstein-Brown formula; the lowest value 1.5, the Meyer-Peter 
formula. For these values, Eq. gives the following equations, respectively 


0.714 
1 


(2) 
n 
m 
Te. 
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For typical contraction ratio B,/B, 1/2, Eq. gives d,/d, 1.64, 
and Eq. gives d,/d, 1.49. The difference these values between 
which view the choice the extreme types sediment transport 
formulas, not very great. 1/4, the difference between 18%-22%. 
Even this difference not much compared with the scatter that common 
the data erosion and sediment transport problems. For the purpose 
constructing theoretical model, the choice sediment transport formula 
will not usually affect the validity and applicability the model seriously, 


OETA. GF THE INLET PROFILE 


FIG. 4.—Plan 2.44-m Long Contraction with Details Transitional Inlet Section 


although advisable adopt the highest maximum value erosion depth, 
e.g., that obtained from Eq. for design purposes. 


Laboratory Channel.—The experiments were conducted rectangular steel 
channel which was 11.4 long, excluding the inlet and outlet sections. The 
width and depth the channel were 0.76 and 0.46 respectively. The 
channel could tilted using hand-operated screw the upstream end. 
tainter gate was mounted the channel about 1.2 upstream the lower 
end the channel order control water depths the channel. Usual devices 
reducing turbulence the upstream end, e.g., honey combs, use smooth 


| INLET CONTRACTION - OUTLET | RK 
| 
FIG. 3.—Plan 1.83-m Long Contraction with Details Transitional Sections 
tele 
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transition the inlet, etc., were used. venturimeter was provided the 
0.20-m diam pipeline for the measurement discharge. The water and bed 
levels were measured using point gages. For the details arrangements for 
measuring longitudinal slope, sand feed, and sand transport, see Ref. 

Bed sizes sand were used the experiments. The average 
size the coarse sand, was 1.53 mm, while the average size the fine 
sand was 0.92 mm. The sands were nearly uniform size; the gradation curves 
are available Ref. 

Long first series experiments, the effective length 
the contraction was 1.83 excluding the upstream (inlet) and downstream 
(outlet) transitions, each 0.46 long (see Fig. 3). The overall length was thus 
2.75 Later, the effective length was increased 2.44 with transitions 
each 0.46 long (see Fig. 4). the 1.83-m long contraction, the upstream 
and downstream transitions were identical shape. The transitional curve was 
simple circular arc, shown Fig. The width the side wall was 
0.127 the total contraction was thus 0.25 channel 0.76 width. 
The angle inclination, was nearly 15°, slightly greater than 12.5° proposed 
Hinds and quoted Smith (4), for smooth transitions. view some 
recent investigations simplified flume inlets (4), 15° inclination was considered 
satisfactory. the very first few experiments, was observed that deep holes 
were scoured near the transition walls the inlet. order check the 
deep scour could avoided sufficiently reduced providing better type 
transition the inlet, the simple circular transition was later changed 
the sinusoidal transition, shown Fig. However, the downstream transition 
was not changed. The local scour near the inlet region was not reduced, and 
continued appreciably deeper than that within the long contraction. 

The contraction the channel was achieved placing wooden side walls, 
each 0.127 wide, within the channel. The walls were placed symmetrically 
opposite each other. The transition pieces, also made wood, were screwed 
the main body the contraction. The whole structure the contraction 
was painted smooth finish. 


Before releasing the water the channel, the sand bed was leveled, both 
the contracted and the uncontracted sections, using wooden screeds 
appropriate dimensions. The sand bed was 22.5 thick. The depth flow 
was adjusted the tainter gate. The water levels recorded the end the 
experiments were used the calculations the flow depths and the bed shear 
stresses, since these values were believed correspond the equilibrium 
conditions bed scour. Each experiment was run for least 

The determination the correct rate sand feed the experiments with 
bed load movement was approximate. Although the rate sediment transport 
was measured the outlet end several experiments, however did not 
indicate the average rate sediment movement the upstream uncontracted 
channel. Visual judgment was used controlling the feed rate during 
experiment. the end experiment, the bed the upstream channel 
was found substantially lowered below raised above its initial flat 
position, then the experiment was repeated with adjusted rate sand feed. 
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The sole criterion steady flow sand was thus the final condition 
the bed the upstream channel. Although very accurate evaluation the 
mean position the bed level the end the experiments cases where 
ripples were formed over the bed, was difficult, fairly good idea about the 
condition sand flow could obtained taking numerous readings the 
bed. Since the sands were uniform size and the ripples were two-dimensional 
and fairly long (0.61 m), the evaluation the mean bed level could 
satisfactorily done. 

The scour depth, d,, shown Fig. was obtained averaging several 
(generally four) depth readings taken along the longitudinal section passing through 
the midwidth the channel, over distance 0.61 The section 0.61 
length was nearly from the inlet and outlet the long contraction. 


(a) 


FIG. Eroded Sand Bed Long Contraction 


For the details sand feed arrangement, measurement feed rate, sand trap 
and measurement sand transport, see Ref. 


Erosion 


Patterns erosion long contraction are influenced the condition 
flow the upstream uncontracted channel. When the bed shear stress 
the upstream channel was less than the critical shear stress, the bed the 
long contraction was eroded, but did not have sand waves over the bed 
the equilibrium condition scour for the bed materials used the experiments. 
the beginning, when the flow was just released, some sand waves might 
appear the long contraction, but their growth was progressively inhibited 
the bed shear stress went continuously decreasing with the development 
scour the long contraction. Finally, the equilibrium condition was reached, 
the sand waves disappeared altogether. 

the case bed load movement the upstream channel, the sand waves 
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may may not disappear the long contraction, depending the flow conditions 
and the bed material size. 


Erosion the long contraction was generally nonuniform, both the 
longitudinal and transverse directions. contour map the eroded bed 
shown Fig. for experiment. Fig. shows 1.83-m long contraction; 
Fig. shows 2.44-m long contraction. Fig. deep local scour near 
the left side wall clearly shown. Two-dimensional sand waves the uncontract- 
upper channel and inside the contraction are clearly shown the figures. 

the inlet region, deep holes were scoured either side close the walls. 
The reason for this was believed the generation secondary currents. 
most the experiments, was observed that there was deeper scour 
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FIG. 6.—Plot Experimental Data /d, and Coordinates for General Erosion 
Inside Contraction 


the region close the right wall the contraction than that the corresponding 
region near the opposite wall. The difference the maximum depth scour 
the two sides some experiments was much 2.5 cm. This was believed 
mainly due higher local velocity near the right wall. The wooden walls 
which the channel was contracted were manufactured accuracy 
0.3 cm. Slight misalinement the contraction walls could have been source 
causing different scour depths the two sides. However, great care was taken 
alining the contraction walls. These walls were removed and replaced 
least three times during the course these experiments, and because the deeper 
scour always occurred near the right side wall, higher velocity near the right 
side wall may have been caused some very slight transverse slope which 
could, however, not detected. Some velocity measurements were actually 
taken which supported the preceding premise. 
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Inside the contraction, the maximum scour generally occurred the center 
the channel close it, although some experiments, the maximum 
erosion was close the right wall. 


When there was bed load movement the upper reach, sand was fed 
appropriate rate that the material removed from the bed surface the 
flow was replaced sand feed. relatively small rates bed load transport, 
the eroded bed the upper reach remained more less flat. When the bed 
load transport was increased, sand waves started forming both the upper 
reach and inside the long contraction. Erosion profiles inside the long contraction 
were similar those when the upper channel sand bed was stable, the difference 
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FIG. 7.—Plot Experimental Data d,/d, and Coordinates for Local Scour 
Near Side Wails Region 


being that the bed was covered with sand waves high rates bed load 
movement. 


Long Contraction.—The experimental data have been appropriately reduced 
order check Straub’s formula for the equilibrium depth scour. The 
results are plotted diagram with and d,/d, coordinates (see Fig. 
6). The prediction from Straub’s formula that the maximum value d,/d, 
occurs 1.0 supported the experimental results. Eq. with 
has been plotted this diagram for the sake comparing the predicted 
values with the experimental ones. Despite obvious scatter, the experimental 
data are well described the calculated values scaled about 20%. Thus, 
the following equation appears give the average value Fig. 
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B, 6/7 B, 1/3 + + -3/7 


The experimental values are fact the values maximum depth scour 
the central plane within the long contraction already mentioned, whereas 
Eq. correlates the ratio the average cross-sectional scour depth and the 
upstream depth with The experimental values are thus expected 
somewhat higher than the theoretical values. Fig. the outer envelope 
curve also plotted, which corresponds the following equation 


B, 6/7 B 1/3 -3/7 
Since the design invariably done against the worst situation, suggested 
that Eq. used design practice. There does not seem any great 


TABLE 1.—Data Scour Depth 1.83-m Long Contraction with Coarse Sand Bed 


“LH and left hole and right hole near the side walls the inlet region, respectively. 


influence the bed material size the scour depth far these experiments 
reveal. Attention must, however, drawn the fact that while the fine sand 
bed reached nearly state equilibrium the coarse bed might continue 
eroding beyond this time limit. Such additional erosion is, however, believed 
relatively small. Since the experiments were performed with only two 
grades sediments, the statement that the sediment size unimportant 
far the maximum relative scour depth concerned, cannot generalized 
unless supported with additional data with other different sizes sediments. 
The average value bed shear stress, for which erosion would occur 
the long contraction, can obtained from Eq. making d,/d, 1.0. 
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Sand 
feed, 
liters grams 
centi- per centi- per 
(1) (3) (4) (5) (6) (7) (8) (9) (10) 
9.60 6.71 15.04 9.14 1.36 1.51 
9.60 8.05 19.86 11.58 1.44 1.26 1.85 
9.60 8.08 22.41 13.11 1.62 1.26 1.82 2.52 
23.30 5.43 15.04 9.45 1.74 0.79 1.30 2.75 
23.30 6.71 19.86 12.50 1.86 0.64 2.27 2.30 2.65 
23.30 7.93 23.83 11.89 1.50 2.27 2.30 2.58 
33.00 3.44 9.93 7.62 2.22 0.90 2.27 2.90 3.58 
33.00 4.27 15.04 8.53 2.00 0.71 3.18 
21.40 3.63 9.08 4.88 1.34 1.26 2.48 3.25 
21.40 4.76 13.34 9.14 1.92 1.00 2.66 3.03 
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The value B,/B, the present experiments was 2/3 that 4/9, 
for which the bed within the long contraction should remain uneroded. This 
value should, however, corrected view the present experimental results 
which indicated that the observed value d,/d, about 1.2 times greater 
than the predicted. After making this correction, 0.37 the maximum 
value which erosion should occur the long contraction with B,/B, 
2/3. the experiments, value 0.63 was reached for which 
d,/d, was nearly 1.38. 

Inlet Region.—There were deep holes eroded the sand bed near the upstream 
transition walls the long contraction either side, already described. 


TABLE 2.—Data Scour Depth 2.44-m Long Contraction with Coarse Sand Bed 


liters 
per 


and left hole and right hole near the side walls the inlet region, respectively. 


The maximum depth scour these holes was nondimensionalized dividing 
the local flow depth d,. Fig. all the data for the coarse, well 
the fine sand, and also the data for the right and left scour holes, are plotted 
d,/d, and coordinates. The following equation also plotted Fig. 
and appears pass through the mean position the plotted points 
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This gives 
2 
Sand 
feed, 
d,, d,, in 
grams 
centi- centi- per 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 
9.60 6.34 15.04 8.53 1.35 1.60 1.60 2.22 
9.60 6.83 16.46 9.75 1.43 1.49 1.75 2.25 
9.60 7.41 19.58 10.67 1.44 1.37 1.85 2.11 
9.60 8.35 21.00 11.89 1.42 1.22 1.94 2.15 
23.30 4.94 15.61 9.45 1.91 0.85 1.09 2.45 2.85 
23.30 5.94 20.71 10.36 1.75 0.71 2.41 3.11 
23.30 6.40 23.55 10.67 1.67 0.66 4.36 2.40 
23.30 4.08 12.20 8.84 2.17 1.03 2.80 3.40 
23.30 6.89 25.25 10.06 1.46 0.61 4.36 2.51 
33.60 3.44 9.93 7.01 2.04 0.85 2.65 3.30 
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TABLE 3.—Data Scour Depth 2.44-m Long Contraction with Fine Sand Bed 


liters 


“LH and left hole and right hole near the side walls the inlet region, respectively. 


order cover the extreme values scour depth, the following outer envelope 
equation also plotted Fig. 


Thus, all the values scour depth actually measured are less than equal 
the calculated values from Eq. 15. The data are given Tables 1-3. 


Straub’s simple one-dimensional model can used for predicting scour depth 
long contraction with reasonable accuracy. For design purposes, the predicted 
depths should, however, scaled about 58% get maximum local 
value. The inlet region close the walls subjected severe local erosion. 
order obtain estimate the extreme maximum depth scour 
this region, the theoretical values should multiplied 2.56. 


Sand 
feed, 
centi- per centi- per 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 
5.20 6.55 12.20 8.23 1.26 1.51 1.45 1.77 
5.20 6.86 16.46 10.06 1.47 1.44 2.05 2.48 
5.20 3.87 7.09 4.88 1.26 2.55 1.58 1.51 
9.04 4.36 9.36 6.71 1.54 1.34 1.82 2.42 
9.04 5.33 12.20 8.23 1.54 1.08 2.18 2.52 
9.04 5.43 13.34 8.84 1.63 1.00 2.20 2.55 
9.04 6.68 16.17 10.67 1.60 0.87 2.06 2.52 
9.04 7.53 21.28 11.89 1.58 0.76 2.25 2.45 
9.04 7.56 21.56 12.19 1.61 0.75 2.27 2.35 
12.30 6.37 9.93 9.45 1.48 1.09 2.05 2.18 
12.30 2.74 5.11 3.96 1.45 1.52 1.55 2.08 
12.30 5.18 15.04 7.93 1.53 0.81 1.09 2.45 2.52 
12.30 4.45 12.20 7.62 1.71 0.94 2.22 2.62 
12.30 5.91 17.31 8.84 1.50 0.71 2.27 2.42 2.58 
12.30 6.55 21.85 12.50 1.91 0.64 4.99 2.18 2.25 
16.60 4.72 14.47 7.93 1.68 0.65 3.18 2.25 2.80 
16.60 6.61 22.70 9.75 1.48 0.47 7.63 2.15 2.45 
52.90 4.05 16.17 6.10 1.51 0.24 20.70 2.25 
62.00 3.90 15.32 5.79 1.49 0.21 27.24 
~ 
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The following symbols are used this paper: 


widths normal and contracted channels; 

size bed material; 

flow depths the normal and contracted channels; 
gravitational acceleration; 

constant; 

Manning’s coefficient bed friction; 

discharge; 

rate sediment transport per unit channel width; 

energy gradient; 

specific gravity; 

flow velocity; 

bed shear stress the normal channel; and 
critical shear stress the bed material. 


g = 
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INTRODUCTION 


the wind driven motion near the free surface carefully examined, one 
often finds sharp velocity gradient that upper few meters water. Van 
Dorn (32) suggests that the drift current depth may little 
one half the surface value. Shulman and Bryson (24) measured very sharp 
gradient and Teeson, al. (28) concluded, based field studies, that the 
velocity gradient near the surface large enough that the draft object 
used determine surface drift will have effect the observation. variety 
laboratory experiments have been performed Shemdin (23), Baines and 
Knapp (1), and Larnaes (14), among others, which indicate large velocity 
gradient exists near the surface, similar logarithmic law the wall type 
profile. 

Fig. based data from Stolzenbach, al. (27) and shows results 
surface drift experiments terms the and the 
The deflection angle the angle the surface drift, measured 
the right the wind speed and the wind factor the surface drift speed, 
normalized the 10-m wind speed. interest note that the 
wind speed common rule thumb for estimating surface oil slick movement. 
Also note that the wind factor the data varies from about indicating 
that the rule oversimplification. 

first try interpreting the data available wind induced currents, 
the writers turned the classical Ekman (9) theory. However, for deep water 
the Ekman (9) theory produces surface current that 45° the right 
the wind direction (in the northern hemisphere). can readily observed 
from Fig. that the measured deflection angles consistently fall between 
and 13° with none close 45°. Several these experiments were examined 
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detail determine bottom unsteady effects could the discrepancy 
between theory and observation but none these were successful. Therefore, 
the most likely explanation would appear the assumption the Ekman 
formulation that the vertical eddy viscosity constant over the depth. 


@i2 
DEFLECTION 


(13) 
(1°) 
TEESON ET AL. (28) 
SMITH 
DOEBLER (8) 


ac 


WIND 
DIRECTION 


FIG. 1.—Comparison Observed Surface Drift Deep Water Ekman Solution 


The vertical eddy viscosity arises from the turbulent frictional force between 
water layers, t,, and t,,. Ekman used the laminar analogy write these stresses 
function the turbulent vertical eddy viscosity coefficient, N,, which 
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HY3 NUMERICAL CIRCULATION MODEL 
can also written terms the turbulent shear stress components or: 


which u’, v’, the turbulent components the velocity the 
and directions, respectively; and the water density. Fig. indicates the 
coordinate system used this paper. 

The component, will zero near solid boundary, and since 
(or will not, general, zero near the boundary, follows that 
must zero the vicinity solid boundary order satisfy the right 
side Eq. waves are not breaking the surface, then the foregoing 
argument implies that must near zero the surface. Thus, based 
physical interpretation N,, concluded that realistic variation would 


shape with near zero the surface and the bottom, and maximum 
somewhere between. 


SEA 
BOTTOM 


FIG. 2.—Definition Variables 


Arguments for vertical variation have been advanced many including 
Fjeldstad (10), Bowden, al. (2), Murray (16), and Thomas (29). Madsen (15) 
assumed the vertical eddy viscosity zero the surface and increase 
linearly with depth. For steady state, his results indicate deflection angle 
approx 10° the surface—a value more consistent with the field observations 
Fig. Other aspects Madsen’s model remain comparatively similar 
Ekman’s model, including the magnitude the surface current and the total 
mass flux. 

Another situation which vertical variation the eddy viscosity results 
more realistic flow pattern, for the case flow long channel 
finite depth with the wind blowing the axial direction. For this situation, 
analytic solution that assumes constant vertical eddy viscosity, yields 
linear velocity profile. instead assumed that varies parabolically (zero 
the surface and the bottom) then the solution for the velocity logarithmic 
profile near the surface and the bottom. This profile substantiated laboratory 
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experiments performed Shemdin (23) and others. 

The model presented here combines advantages both analytical and 
numerical approaches. particular, the model: (1) Makes relatively few simplify- 
ing assumptions; (2) has reasonable computation requirements; and (3) yields 
continuous vertical variation for the velocities. The model logical extension 
model presented Heaps (12) which uses the Galerkin weighted residual 
technique simplify the Navier-Stokes equation. However, the model presented 
this paper does not assume the vertical eddy viscosity constant 
the vertical direction. its present form, the model being successfully applied 
several estuarine and coastal areas for the calculation wind and tidally 
generated currents. 


many sea and lake circulation problems practical importance, the vertical 
velocity component and its gradient are considerably smaller than the horizontal 
velocities and gradients, and can neglected. With this assumption, the 
momentum equation reduces expression the hydrostatic pressure 
distribution. The component the momentum equations used the model, 
with variables defined Fig. is: 


water depth some datum (in this case, low water); gravity; the 
Coriolis parameter sin where latitude; and angular velocity 
the earth; atmospheric pressure; water density; subscript designates 
the surface value; and horizontal eddy viscosity. similar expression 
used for the direction. The vertically integrated form the continuity 
equation for incompressible fluid 


The surface boundary condition is: 


which and the surface shear stresses induced the wind the 
direction and direction, respectively. 
The shear stresses the bottom boundary can expressed 
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which and denote the components bottom friction the and 
directions, respectively. 
linear bottom friction assumed, 


The governing equations are simplified using expanded version 
Galerkin weighted residual approach used Heaps (11, 12). detailed 
description the development will not included here but can found 
Cooper and Pearce (4) Pearce, al. (20). However, since the development 
somewhat unusual, appropriate include brief outline. 

The application weighted residual technique begins assuming functional 
relationship approximate the true solution: 

d,Q, 

I=1 
which and the “‘trial functions”; prescribed function; and 
coefficients the prescribed functions which are determined the solution 
process; the number terms needed reach the required degree 
accuracy; and and terms which are functions but constant 
time. 

Heaps (11,12) chose cosines for the prescribed functions. Other functions 
have been tried Davies (7). Cosines have the advantages being well behaved 
when differentiated and forming orthogonal set. 

Using cosines, the final form the trial functions was chosen 


p,H N, p,a N, i=l H 


2 
a,z 


which all symbols are specified constants constant functions except for 
the undetermined parameters, and d,. The variable has the value 
the bottom. The values are determined satisfying the bottom boundary 
conditions and can found the solution to: 


and are the results evolutionary process the writers. Originally, 
simpler form (4) without the logarithmic term was chosen. was discovered, 
however, that using variation which became small near boundary, 
convergence problem occurred the model. was reasoned that since the 
velocity profile the upper several meters water develops rapidly, then 
form for u,, that contains logarithmic term would aid convergence. The 
and shown and dramatically improved convergence that, 
general, only three cosine terms the series are needed for the variety 
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applications investigated far. no-slip boundary condition the bottom 
was included the initial formulation and caused similar problem with 
convergence. However, the solution this case could not remedied 
the case the surface and slip bottom boundary was imposed, Eq. 

The trial functions, Eqs. and are substituted into the momentum equations, 
Eq. and its analogy. general, there will error associated with 
this substitution, since the trial functions are not the exact solutions the 
momentum equations. This error sometimes referred the 
and minimized specifying that the integral the residual over region 
zero. The Galerkin method consists multiplying the residual and 
integrating over the water column. The integration yields: 


which the term has been differentiated and conveniently separated. There 
equivalent Galerkin expression for the direction. 


FIG. 3.—Functional Form for Vertical Variation 


Observe that the lower limit the integration Eq. has been moved 
would become function time, thus, complicating the solution and increasing 
computation time. 

order evaluate the integrals, the vertical variation must specified. 
This done approximating with series linear segments shown 
Fig. The slope a,, the intercept B,, and the number segments are 
assumed given. [When more than one linear segment used (i.e., then 

Carrying out the integrations implied Eq. 11, transforms the component 


Eq. into equations terms the undetermined parameters and 
or: 
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—-eE C,;- + + c or 9 


which A,, B,, and E,, constants which arise from the integration; and 
the terms considered Neumann (17) and Neumann and Pierson 
(18). equation similar Eq. exists for the direction. the present 
formulation, the convective terms (those the form and 
density gradient and atmospheric pressure gradient terms have been neglected. 

Note that the constants B,, A,, and E,,, need evaluated only once 
when changes time. The time scale for significant changes 
the order hours, whereas the time step the order minutes. Thus, 
the computational cost the evaluation the constants relatively small. 
Addition the vertical eddy viscosity variation adds the E,, term. Without 
this term, the number calculations per time step per grid element would 
increase linearly with (recall the number cosine terms the series 
expansion, see Eq. 7). With the term, the number calculations will increase 
approximately squared. 

Eq. and its equivalent the direction represent set equations 
with unknowns (c,,d,,n). solve for the unknowns, one more equation 
linking c,, d,, and must specified. This accomplished using the continuity 
equation, Eq. Expressions for and Eq. can obtained integrating 
the trial functions over the depth. The variables, and are substituted 
into Eq. which will discretized along with Eq. and its counterpart. 
Thus set equations now available deal with the set 2/’ 
unknowns. 

The method chosen for the horizontal discretization the finite 
difference scheme used many modelers such Platzman (21) and Reid 
and Bodine (22). Originally selected for simplicity and vehicle for examining 
the suitability using the Galerkin method, the split-time method has performed 
well but should not considered limiting. 

Numerical three-dimensional (3-D) models can limited scope because 
their computer time and storage requirements are too large for problems 
practical interest. One the major objectives this research was develop 
three-dimensional model which could applied inexpensively. The number 
multiply-divides carried out the program per time step per grid element 
programmed, explicit, vertically-averaged circulation model. Weare (33), esti- 
mates about 200 multiplications per grid point per time step for one vertically- 
averaged, explicit, finite difference scheme. 

The split-time finite difference scheme conditionally stable. expression 
indicating the stability criteria for the full model described here has not yet 
been derived. Two conditions serve provide reasonably good choice for 
the time step, At. The condition [Courant, al. (5); Platzman (21)]: 


which the horizontal grid spacing, must satisfied. general, Eq. 


HY3 
(13) 
2gH 
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the only condition needed for vertically-averaged explicit models. has 
been found, however, that Eq. alone not sufficient for this model; 
additional condition needed. approximate form can obtained examining 
the Galerkin statement, Eq. 11. Ignoring the horizontal eddy viscosity terms 
and the Coriolis term, assuming constant, and carrying out the indicated 
integration produces: 


Performing forward difference time, with time step yields: 


The method von Neuman [Carnahan, al. (3)] introduces c,(t) 
the harmonic component perturbation. That substitution yields upon 
rearrangement: 


which the amplification factor, which must obey the criterion, 
yield stable scheme. Since the second term the right side Eq. 
always positive, the limiting condition for is: 


ANALYTIC 


a 
= 
= 
= 
= 
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FIG. 4.—Comparison Between Model and Analytic Solution for Case Finite Channel 
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Experience operating the model far indicates that Eq. too lax criterion 
and that 


gives better estimate the time step. The discrepancy undoubtedly due 
the many simplifying assumptions made during Eq. 17. While Eq. provides 
working criteria, hoped that more complete analysis can developed 
the future. 


Test Cases 


address questions concerning convergence, stability, and verification 
computer coding, comparisons were made analytical solutions. The problems 
for which closed form solutions can found tend for simplified geometry 
and flow situations. Nevertheless, this aspect the process valuable since 
provides the only situation where nearly perfect agreement can expected. 

the case wind over closed-ended channel with linearly varying eddy 
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VELOCITY (M/SEC) 


FIG. Between Model and Analytic Solution for Case Linearly 
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FIG. 6.—Comparison Between Model and Shallow Water Ekman Solution 0.094, 


c, = 


viscosity, flow the surface along the direction the wind expected 
with return flow along the bottom. The solution given Pearce, al. (20) 
and the model results are compared with the analytic solution Fig. The 
comparison good with only three cosine functions (i.e., 3). illustrate 
convergence, the solution for five cosine functions also shown Fig. 
Here and elsewhere, three cosine functions appear provide sufficient accuracy. 
Many one-dimensional (1-D) comparisons the type shown above have 
been carried out and have shown similar characteristics rapid convergence. 
Comparisons 3-D cases have also been made. Fig. shows comparison 
between the numerical model and the analytic solution for the case wind 
induced flow water body constant depth and infinite lateral extent. 
The problem basically the same the steady-state Ekman solution (9) except 
that this case assumed vary linearly. The figure indicates essentially 
identical comparison using only three cosine terms. Similarly good results are 
obtained assumed constant the vertical indicated Fig. 


The comparisons described the previous section confirms the model 
formulation and coding for simple cases, but further investigate the charac- 
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teristics the model and gain insight into the functional form the vertical 
eddy viscosity parameter, series comparisons laboratory data was 
undertaken. Experiments suitable for comparison have been carried out 
Shemdin (23) and Baines and Knapp (1). Before making comparisons this 
data, however, some variation for must decided upon. Furthermore, 
select method algorithm for choosing values for and c,. 

Distribution distribution with depth is, course, matter 
considerable interest and has been addressed variety ways. Based 
Refs. and others, reasonable distribution for would appear 
shown Fig. The justification for that distribution described 
the following paragraphs. 


fe -- 


FIG. 7.—Distribution Used Model Simulations 


Bowden, al. (2) inferred the variation from tidal current observations. 
Although the data showed considerable scatter, concluded that the values 
tended highest the middle the water column and lowest top 
and bottom. This result consistent with the earlier argument the Introduction. 
Bowden’s recommendation use model based constant eddy viscosity 
the bulk the fluid above friction layer near the bottom. cases where 
wind generated currents are important, boundary layer near the surface 
also important and thus, from Bowden’s recommendation, friction layer near 
the surface also justified. 

Csanady (6) provides formula for the eddy viscosity boundary turbulent 
interface and covers the law the wall region near the interface where the 
velocity. distribution logarithmic. The eddy viscosity constant the 
region very near the interface. The core region equivalent the laminar 
sublayer near solid boundary except that may present 
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which increase the effective eddy viscosity value greater than the molecular 
value, Fig. indicates the various regimes. Csanady found the depth 
the core layer generally the order for one experiment, 

The problem remains find outside the core region. Townsend (31) has 
compiled data the subject, including turbulent flow pipes and channels, 
and notes that the simple assumption that the eddy viscosity constant 
outside the wall layer describes the velocity distribution for various flows. From 
similarity arguments, Townsend concluded that the wall layer 


which the distance above the bed; and von Karman’s constant, 
0.4. This expression satisfies the law the wall. 


Core Layer 


Wall Layer 


Constant N 
Vv 


FIG. 8.—Various Regimes 


For the case open channel flow, Townsend indicates that the limit 
the wall layer occurs the position where: 


which the channel depth; and the flow Reynolds number 
Other sources indicate values about the same magnitude. The flow Reynolds 
number, 12, chosen here for convenience. assumed that 
much smaller than the constant region, then Eqs. and 
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can used evaluate the distance, z,, from the surface which becomes 
constant or: 


Taking 0.4 and the flow Reynolds number, implies that 0.2H. 
Thus, for open channel, the wall layer comprises about 20% the water 
column. Following this analogy for the surface, provides distribution 

Choice c,.—The boundary layer the bottom simulated the use 
slip coefficient, c,, the bottom. The bottom stress given the numerical 
model pc,u,, which the velocity the bottom. this equated 
the shear stress derived from quadratic formula such the Darcy equation 
for can obtained reasonable value can estimated. 


Comparisons Data 


Baines and Knapp.—The experiments Baines and Knapp (1) were performed 
air motion produced shear stress the surface over the flume. 
Two air flow conditions were used with average velocities 6.1 m/s and 
3.9 m/s, respectively. For the 6.1 m/s case: 2,870 which 


depth 0.328 m). The corresponding value 0.00875 m/s and 
from Csanady (6) the effective core surface eddy viscosity, 2.5 
The eddy viscosity the bottom the wall region then 
for that experimental setup. The solution not sensitive c,. Fig. 


Data; 
Model; Theoretical 
Model; Best Fit 


FIG. Between Model and Laboratory Experiments Baines and 
Knapp 2,870) 
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comparison the model output the Baines and Knapp data. The parameters 
chosen above provided the theoretical solution. Except the return flow 
jet near the bottom, the fit quite good. This behavior expected 
since the bottom boundary layer not being modeled. better fit the surface 
can seen with the fit’’ line Fig. The only change that 
has been lowered m’/s. More than one possible explanation for 
the difference the theoretical and the are apparent. Small 
errors could produce discrepancy. The estimates from Csanady (6) 
are approximate. Finally, the presence waves the surface affects N,. 
Slight changes the experimental procedure, such the addition detergent 
Baines and Knapp, may affect the wave climate the apparatus and, therefore, 
N,. 

The fit the 1,900 case Fig. similar. From Csanady, obtain 
2.0 m’/s and this compares the best fit case 1.0 


1900 


Model; Theoretical 


Model; Best Fit 


FIG. 10.—Comparisons Between Model and Laboratory Experiments Baines and 
Knapp 1900) 


m’/s. The analysis based the flow Reynolds number provides 1.5 
m’/s for the interior the fluid which compares well the data. 
about 0.03. 

Shemdin (23) were performed laboratory flume 
deep, long, and 1.83 wide. The flume was divided into two channels 
equal width. One channel the flume had wind tunnel over induce 
water motion. The other channel was used provide path for return flow. 
beach covered with stainless steel shavings was located the end the 
wind tunnel. Three distinct runs were made with wind speeds 9.1 m/s, 
5.4 m/s, and 3.1 m/s. The input for the numerical model was designed 
that the surface shear stress agreed with the shear stress measured Shemdin. 

order model the situation Shemdin created his experiments, boundary 
condition was established the model which allowed for net flux the 
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upstream and downstream boundaries. The flux the amount water being 
driven through the stainless steel shavings and along the return channel and 
was determined from the Shemdin data. 

For the case shown Fig. 11, the theoretical model gave 
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FIG. Between and Laboratory Experiments Shemdin 
9.1 m/s) 


Vv 
= 
a 
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FIG. 12.—Comparison Between Model and Laboratory Experiments Shemdin 
5.4 m/s) 


was chosen 0.02 m/s. The comparison very good, differing the surface 
about cm/s. 

The comparison for the 5.4-m/s case similar. The theoretical model 
0.01 m/s. Fig. shows that comparison. The velocity the surface 
matched very closely. The return flow indicated about Fig. 
not modeled. difficult speculate about the apparent three-layered 
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nature the laboratory data, although natural thermal statification seems 
distinct possibility. comparison the case 3.1 m/s led comparable 
results near the surface, but the data indicated even stronger three-layered 
nature. 

The values chosen for correspond unreasonably high values Lower 
values for corresponding more reasonable values were used, but 
the numerical model continued oscillate for many hours. conjectured 
that the steel shavings the end the channel and the return channel effectively 
increased the friction the system. These effects were included increasing 
the values indicated previously. 


Based compilation current data the surface, concluded that 
the classical Ekman theory not appropriate for calculating wind generated 
currents near the surface. The cause found the assumption Ekman 
vertical eddy viscosity, which constant over depth. 

numerical model described briefly, which allows for variation 
over depth, well the and directions. The model was developed 
using Galerkin technique first introduced Heaps (11). Heaps’ method 
expanded that the model can include sharp, wind driven profile near the 
surface. The method computationally reasonable, requiring about the same 
number computations vertically averaged formulation. Simple stability 
criteria are presented. 

order test coding and accuracy, the model compared several analytic 
test cases. relationship for choosing versus depth presented. That 
relationship based flow Reynold’s number u,H/N,, from Townsend 
(31) and core value from Csanady (6). Choice the bottom friction 
coefficient, c,, based from the Darcy equation. 

Experimental data Baines and Knapp (1) and Shemdin (23) are compared 
the numerical model results. The variables and are shown the 
method described previously. The numerical model provides good comparisons 
the measured data. The writers are currently making comparisons field 
data the Gulf Mexico and other areas with encouraging results. 
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The following symbols are used this paper: 


bottom friction coefficient; 
undetermined coefficients; 
Darcy’s also Coriolis parameter; 
gravity; 
water depth; 
number functions used summation; 
horizontal and vertical eddy viscosity; 
atmospheric pressure; 
Reynolds number; 
subscripts designating surface, bottom; 
and components water velocity; 
turbulent velocity components; 
friction velocity; 
and flux components; 
Cartesian coordinates; 
distance from surface region constant 
coefficients for distribution; 
time step; 
water depth above datum; 
Von Karman’s constant; 
water density; 
shear stress, surface shear stress; and 
trial function. 
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MILLENNIAL CELEBRATION KARAJI’S HYDROLOGY 


One thousand years ago, the brilliance Persian scientific and cultural activity 
was dominant. While scholarly production was extensive, with many contributors, 
attention herein drawn the writing and transcription the unique scholar, 
Mohammed Karaji. celebrate and acknowledge Karaji honorable, especially 
the light the recent claim revealing the discovery his original manuscript 
hydrology. 

This paper presents brief account Karaji’s contributions science, 
especially ground-water hydrology. The correlation between Karaji’s thought 
and knowledge with later day scientists brought forth. Response the 
question—Was Karaji ahead his well-illustrated throughout the 


article. Because the discovered manuscript and its assumed original nature, 
the view the writers that the material should not neglected lost 
among the plethora technical matter, but should duly noted and brought 
the attention science historians and hydrologists. 


Genius 


Mohammed Karaji was Persian scholar the late tenth century (12). While 
identified also names Karadji, Karagi, Al-Karaji, Hassen-al-Hasib, 
his complete name was Muhammed ibn al-Hasen al-Hasib al-Karaji (16,18). 
the mathematical world, Karaji was more often referred Al-Karkhi 
(3,20). Regrettably, much the production, brilliance, and contribution Persian 
scholars during their period great scientific flourishment, from the ninth 
century through the twelfth century, has been lost due language barriers 
and errors the etymological description (as the stresses inflection within 
Arabic letters) translated works (3,20,22). 
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Muhammad ibn Musa al-Khwarizmi, called Alkhwarizmi, wrote the first 
comprehensive book the field algebra (c. 825 A.D.). The book, written 
the Arabic language, was titled Hisab al-jabr (and was referred 
simply Al-Jabr). Latinized form, the title reads Liber algebrae 
almucabola, i.e., The Science Reduction and Confrontation. Note that the 
Latinization the Al-Jabr name led the word algebra; also, the Latinization 
the Alkhwarizmi name gave the word algorithm (3,10,11). Although Alkhwariz- 
mi’s work was the first Arabic book the subject algebra, was presented 
without symbolism. Karaji wrote superior algebra text (in Arabic, the language 
the times), the first with symbolism (5) (c. 1000 A.D.), which provides the 
basic foundation and principles known today. Karaji’s text was titled Fakhri 
honor his patron and protector, the vizir (minister) Fakhr al-Mulk (10,12). 
work resents the laws exponentiation, the operations with surds, e.g., 

addition law radicals, e.g., also novel 
method (apparently method known today finite mathematics induction) 
summing the cubes the natural numbers, i.e., showing that 
was known Nicomachus (c. 100 A.D.); also generalized and solved Diophan- 
tine-type problems ingenious use parameters (such problems determine 
cubic equations through geometrical means, etc. (20,21). 


Transmission Genius ENGINEERING 


engineering, Karaji personally invented ingenious and exceptionally exacting 
instruments that proved useful the surveying and the tunneling ganats 
(underground channels for conveying water). Perhaps most importantly, Karaji 
expounded the basic principles hydrology. Evidently Karaji was familiar 
with the principal gravitational attraction, although did not formally state 
analytical formulation [as today’s useful Newtonian law (14)]. Further, 
appears that Karaji understood the proportionality law between the force 
gravity and the mass body. Such descriptions, well extensive details 
surveying and the hydrology ground water, are found recently discovered 
manuscript Karaji, entitled (in Arabic) al-Miyah al-Khafiya, which 
translates The Extraction Hidden Waters the Surface (12,18). 1966, 
Persian translation this manuscript was produced; today, writing projects 
are underway translating said works into several additional languages 

authenticated and accepted historians science, Karaji’s manuscript 
would serve the oldest text the subject ground-water hydrology. 
the search for (ground water), the manuscript delineates Karaji’s 
universality with physical laws, illustrated the quotation, groundwater 
moves through cracks and cavities from points farther from the earth’s center 
points closer it, and that every object has this tendency, with the property 
becoming stronger for heavier objects (translated from transcription 
the manuscript) (12). 


History 


The origin ground water has been the subject speculation from the 
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earliest days documented history. During the little attention 
was devoted the subject. Many scholars followed Biblical thoughts concerning 
the origin water. Philosophers such Aristotle and Seneca believed that 
springs could not derived from rainwater alone. the and 
early years the Renaissance, scientists such Kepler, Kircher, Magnus, 
and Schottus exhibited great interest the subject and expressed their notions 
form that essentially different than the basic knowledge today. Kircher, 
example, with great appetite for imagination, described the earth (in 1619) 
medium which, drinks water from the sea and digests like 
beast and, result this metabolism, springs Kircher further 
attributed the mystical whirlpool the the coast Norway, 
the mouth underground caverns (1,2,9). 

Thus must appear that explanations offered early philosophers the 
origin ground water were possibly erroneous. Perhaps strange, but 
explanatory notes for some these beliefs are attributable to: (1) The desire 
God; (2) the theory that the surface the oceans are higher than the surface 
land; and (3) the belief that rock pressure water upward through 
underground caverns, etc. Leonardo Vinci, outstanding figure the 
Renaissance and genius the western world, was one the earliest present 
true explanation the origin rivers. Bernard Palissy, writer the 
sixteenth century, appears the first recognize that rain and snowmelt 
serve the source springs and rivers. Later, 1715, Antonio Vallisnieri 
emphasized that rivers and springs originate their waters from rain, and also 
described the mechanism free-flow operations artesian wells, presenting 
extensive graphics well cross sections (1,7,9). 

Meinzer states, prior the latter part the seventeenth century, 
was generally assumed that water discharged from springs could not derived 
from rain’’ (15). was also generally believed that rain was not adequate 
quantity, and that the earth was too impervious permit the penetration 
rainwater far below the surface earth. Two the main hypotheses extended 
Greek philosophers persisting with confusion through the middle the 
seventeenth century were: (1) That sea water was connected springs via 
subterranean channels; and (2) that air was condensed into moisture cold 
caverns inside the earth. was not until 1686 that scholars like Mariotte introduced 
correct ideas regarding the source springs and the infiltration rainwater 
through the ground (7,9). Much later, 1921, Ototsky, Russian hydrologist 
rejected this infiltration concept reasoning that deep infiltration occurs only 
exceptional circumstances within special regions, and further assumed 
that condensation serves important factor the origin ground water 


(15). 


During the Middle Ages, little recognition was afforded scientific writers from 
nonwestern nations because unfamiliarity language and inaccessability 
manuscripts (8,16,17,24). Many physical concepts that have been described 
basic fundamental learning western scholars the seventeenth century 
and eighteenth century were apparently known Persian scholars centuries 
earlier. Karaji, cited previously, knew the origin water,”’ i.e., 
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ground water, its occurrence, and its relation rock formations, aquifers, 
and the hydrological cycle. Karaji also recalled that few dozen deep-rooted 
plants (phreatophytes) serve excellent indicators the occurrence ground 
water. Biruni (973 A.D.-1048 A.D.), another masterful and brilliant Persian 
scientist and philosopher referred whose contribution 
amounts some 180 outstanding publications various fields (19), was also 
among the first give correct explanation the mechanism artesian 
wells (7,22). Such work, however, remained unnoticed and was repeated seven 
centuries later Vallisnieri (9). 

Karaji, contemporary Biruni, demonstrated interest the practice 
engineering addition his promotion science, and displayed his talents 
and knowledge through the subject The following excerpts 
from the recently discovered manuscript Karaji present brief indication 
his contributions groundwater science (12). These excerpts are translations 
the authors, from Persian English. 

Karaji states that 


the source springs and all water bodies inside the earth, well 
all surface waters, comes from rain and snow the rains 


and snows cease, then waters are depleted and the earth gradually 
wasted. 


Also 


Rains and snows percolate the earth with excess waters running off into 
the seas, and that percolation ceases when hard clays other solid 
obstructions are encountered. 


Such vivid delineations the existence groundwaters point the description 
processes which today defines the hydrological cycle. 
Further, Karaji states (12) 


heard that large springs fresh water exist some islands. 
The source origin these undoubtedly not the water from the seas 
that surround the islands, because the water surface the seas lower 
than the surface the islands; rather, this source located farthest 
points and higher elevations. 


the relation geological formations groundwater occurrence, Karaji 
classifies rock formations their color—black dark rock mountains typifying 
the existence the richest groundwaters, also, that the abundance ground 
water decreases, respectively, for green, yellow, and red rock mountain forma- 
tions. This subject was totally ignored until valuable contribution was made 
the nineteenth century Gabriel Daubree, French geologist, the relation 
geological structures the occurrence and movement ground waters (15). 
Karaji noted that ground waters, analogous surface waters, can stagnant; 
also, that still waters can located under vast deserts and lowlands. Thus 
Karaji was illustrating his knowledge ground-water movement. 

Further evidence the extent Karaji’s knowledge ground-water movement 


HY3 HYDROLOGY 
exhibited the following examples (from his manuscript): 


The water levels wells community near the Diayla River rise whenever 
there increase river flow, and drop with decrease flows 
such that the water level the wells the same that the river. 
With levees (or dikes) constructed along the river, prevent flooding 
adjacent lands, the water levels the river rise higher and produce 
corresponding rises the level water wells, consequently flooding 
homes the community. The water from the river was fresh, though 
muddy, while clear but saline the wells. 


Karaji thus reasoned that this flow was due the absence obstructions 


barriers, i.e., impermeable layers earth, stop the flow 
Further, Karaji wrote 


From rock, projecting out river near the village Saveh Province, 
water spouts through vents. The river water fresh, but the spouted 
water bitter that induces ill effects drinkers. The source 
this bitter water certainly not the river, but the percolating waters some 
distance away which its taste traveled through soils enroute 
the rock. 


These examples indicate Karaji’s knowledge unconfined and confined 
flows, respectively. 

Water quality also considered Karaji’s manuscript. Karaji presents methods 
determine such water qualities hardness, turbidity, taste, color, and odor, 
also presenting ingenious means which water quality can improved (as 
lowering hardness and reduction brine content). show his 
the physical laws which relate pressures and volumes gases, the process 
extracting fresh water from the bottom shallow seas sailors detailed. 
The tools required for such extraction include heavy lead vase with small 
holes the bottom and slim neck the top connected long thin tube 
fine waxed leather, small globe fitted snugly the bottom the thin 
tube, and two long cords—one extending from the sea vessel the handle 
the vase and the other running through the tube and fastened the inserted 
globe. extract fresh water, the entire apparatus dropped into the sea 
and allowed touch the bottom. pull the cord attached the globe 
results increase the volume air creating suction (thereby reducing 
the pressure the vase) and allowing fresh water from the ground beneath 
the sea seep through the small holes the bottom the vase. The container 
then pulled the surface and the routine repeated until sufficient 
quantity fresh water extracted. 

the subject qanats, underground aqueducts, Karaji demonstrated 
his superior technical talents. While some geohydrologists have referred the 
construction most extraordinary works ancient man’’ (23), 
the methods construction described Karaji are essentially the 
same those practiced today centuries later. The fact that tunneling 
requires extensive skills surveying, the provision adequate ventilation, 
and the prevention mining disasters shows the extent technological 
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advancement the Persians qanat construction (6). Indeed, the planning 
and establishment many villages, cities, and centers learning followed 
such construction technology. Tragically, these centers were sacked and destroyed 
invading warriors. Consequently, much the treasured scientific knowledge 
was destroyed lost. 


The history scientific endeavors reveals that Karaji was productive and 
highly knowledgable areas mathematics, ground-water hydrology, surveying, 
and the engineering tunneling. Karaji may perhaps have shared many important 
scientific ideas with his brilliant contemporaries (namely, Biruni and Avicenna). 
However, errors ethnography and unfamiliarity with inaccessibility other 
Karaji works other than those mathematics have caused the omission 
Karaji’s name applied scientist western writers. The achievements 
Persian scholars, which flourished during the ninth century, the tenth century, 
and the early part the eleventh century has indeed formed brilliant era 
the history scientific learning. 

Continual political turmoil and religious factionalism precipitously led periods 
inconstant scientific production and lack high degree uniformity 
learning. From the middle the eleventh century, science efforts Persians 
and Arabs were reduced such state decline that never again rose 
the scholarly level attained Karaji. Overrun invasions during the twelfth 
century, centers learning were sacked and manuscripts were destroyed. 

Additionally, western scholars have praised Persian scientists admirable 
custodians preservers world intellect, which transferred into Europe and 
the western world. However, western scholars have questioned the originality 
Persian writings. Persian scholars made extensive translations Greek and 
Roman classical works (such the works Apollonius, Diophantus, Euclid, 
Vitruvius, etc.). 

Regardless the stature assigned Persian scientific endeavors, Karaji stands 
alone prominent contributor the hydrological sciences. The recent discovery 
Karaji’s manuscript identifies the manuscript the oldest 
known text ground-water hydrology. The reception and authentication 
this manuscript science historians should strengthen the value this 
millennium work. Writers history science should thus include the significant 
Persian contributions during the tenth century and eleventh century. also 
very appropriate call Karaji founder this millennial 
celebration. 

conclusion, Karaji viewed the subject ground water with such great 
importance that the preface his manuscript reads: 


not know any profession more useful than that the extraction 
ground water the surface. Because this profession, land reclaims 
and the liveliness man finds organization and plenty benefits. 
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TURBULENCE MEASUREMENT STUDY 
Franklyn Griffith, and Charles ASCE 


Increasing concern regarding the environmental effects changes flow 
patterns for bodies water and streams demands that flow parameters charac- 
terizing turbulence measured and understood. necessary, therefore, 
determine the significant parameters and devise measurement methods that 
utilize instrumentation compatible with field use. Discharge and current are 
generally used describe flows, with dissolved oxygen, sediment load, tempera- 
ture, etc. used characterize the content condition flow. Rapid mixing, 
sudden pressure changes, and sharp pressure gradients that are associated with 
turbulent flow have impact sediment transport, clustering soft organisms, 
and channel stability. 

This paper presents the results tests intended evaluate the feasibility 
using modified commercially available electromagnetic current meter 
measure parameters which would characterize turbulent flows from field mea- 
surements rivers and lake inlets. anticipated that the techniques described 
here will used the analysis the flows channels connecting Lake 
Ponchartrain and the Gulf Mexico. 


electromagnetic flow meter, Model 511 manufactured Marsh-McBirney 
Inc., Rockville, Md., with modified time constant circuit was used basic 
sensing instrument these experiments. Data were recorded Gould (Brush 
222) two-channel oscillograph recorder and Triplet Model 603 meters were used 
for voltage measurements. Oscilloscope photographs were also used laboratory 
evaluation. Turbulent outflow the Inner Harbor Navigation Canal (IHNC) 
Lock, New Orleans, and tidal flow under the Seabrook Bridge New Orleans 
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and flow the Mississippi River New Orleans were used field test situations. 

The flows are characterized the following parameters: (1) 
average root mean square (rms) the fluctuating (2) time average 
velocity resultant (D) averaged over several seconds; (3) spectral composition 
from Fourier analysis strip chart recordings; and (4) the 
the modified turbulence intensity. 

Previous methods determining the foregoing parameters present problems 
when attempts are made apply them the field where number crews 
must use the instruments from small boats, bridges, and locations where access 
difficult. The instrument package described here relatively light, battery 
powered and rugged. Techniques described the literature such the laser 
doppler anemometer (1) are not portable and require delicate alignment. The 
hot film sensors used successfully McQuivey (3,4) are delicate, etched 
abrasives carried the sediment load and easily contaminated. 


FIG. Probe for 511 Marsh-McBirney Electromagnetic Current 
Meter 


The Model 511 was evaluated potential candidate for yielding turbulence 
information. The probes for this instrument are shown Fig. They consist 
electromagnet which modulated Hz, two pair sensing electrodes 
orthogonally mounted 1-1/2 in. (3.8 cm) diam ball. Circuits and physical 
details are contained the instrument manual. These systems have been used 
successfully field discharge measurements and fixed flow monitors 
remote locations for several years and have proven reliable. The instrument 
has selectable velocity scales fps (3.05 m/s) and time constants 
sec with the two orthogonal velocity components simultaneously displayed 
panel meters. separate output with representing (0-m/ 
flow velocity available. This output was recorded during the 
experiments described here. 


Modification the output circuit consisted reducing the time constant 
replacing the capacitor (C12) with capacitor. This increased 
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the frequency response approx Hz. The instrument’s frequency response 
shown Fig. was determined the use oscilloscope measure 


RELATIVE AMPLITUDE 


FREQUENCY 


FIG. Response over Range Interest 


FIG. 3.—Example Output Waveform for 1-Hz Simulated Velocity 


the peak peak output voltages and also high impedance a-c voltmeter 
indicate rms values. Results from both methods were similar and showed the 
3-dB point Hz. Beats between the modulation and the magnet power 
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frequency were evident near but were lower amplitude and did not 
contribute significantly the results. 

Figs. and show examples instrument outputs for simulated velocity 
fluctuations and Hz, respectively. This signal form was observed 
for frequencies approx Hz. Recordings such the ones shown were 
made inducing electric potential into the water containing the probe. 
This potential consisting the amplitude modulated 30-Hz carrier which was 


FIG. Waveform with Ground Loop Noise Shown Lower Trace; 
Upper Trace Shows Potential Waveform Used Simulate Fluctuating Velocity 


phase locked the magnet driver and simulated the signal expected from 
sinusoidal velocity fluctuation from turbulent flow. necessary look 
the effect the instrument and the analysis methods the spectral distribution 
measurements. 

Oscillograph recordings were made frequency increments from Hz-14 
for various noise conditions. number these were digitized 0.03-sec 
increments and Fourier analyzed determine the form output voltage for 
various frequencies simulated sinusoidual velocity fluctuations. Two methods 


FIG. 4.—Example Output Waveform for 4-Hz Simulated Velocity Fluctuation 
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digitizing these recordings were tested determine the type transfer 
function and noise effect for each. The envelope for high noise condition 
shown the lower trace Fig. Use the short-term averaging the 
signal for digitizing removed the 30-Hz component and yielded sharp lines 
the modulation frequency. The Fourier analysis this average for 
Table shows sharp line with significant spurious components. Information 
therefore best obtained Fourier analysis the short-term (over period 
about 0.05 sec) average. 

The noise contribution difficult assess since depends measurements 
made flowing but nonturbulent medium. There always uncertainty unless 
special care taken prepare and establish laminar flow. Upper limits, however, 


TABLE 1.—Fourier Components Output with 4-Hz Simulated Turbu- 
lence Fluctuations Induced Probe Tank 


Frequency, Hertz 
(3) 


FIG. Longitudinal Component Turbulence Recorded Laboratory 
Test Tank 


can obtained using the best steady flow available and comparing with 
values calculated methods described the manufacturer’s manual. The 
wide-band electronic noise given rms fps, which the output 
time constant, seconds. The system, modified, would yield broad-band 
noise 0.1 rms fps (0.03 m/s). This noise water flow which appeared 
steady approx fps (0.305 m/s) was measured 0.13 rms fps (0.04 
m/s). evident from recordings and oscilloscope observations, that the 
short-term averaging will greatly reduce the noise components near Hz. The 
noise interference the digitized data will thus reduced considerably below 
the levels indicated which are quoted for broad-band measurements. Since the 
levels fluctuating signals turbulent flow can high five times greater 


(1) (2) 
2.74 0.00 
0.06 0.11 
0.05 

—0.06 0.01 

—0.08 0.26 
0.02 —0.08 
0.03 
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TABLE 2.—Fourier Components Obtained from Longitudinal Turbulence Component 
Generated Test Tank 


Frequency, 
Hertz 


FREQUENCY 


FIG. 7.—Frequency Spectrum Longitudinal Turbulence Component Observed for 
Turbulent Mixing Test Tank 


(1) (2) (3) 
6.71 0.00 0.0 
1.25 0.5 
0.50 0.35 1.0 
0.42 1.5 
0.29 0.49 2.0 
0.37 0.42 2.5 
0.12 0.94 3.0 
3.5 
4.0 
0.15 4.5 
0.21 5.0 
0.29 5.5 
0.02 0.08 6.0 
0.12 6.5 
0.05 7.0 
0.14 
4 
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than the broad-band noise, the measurement turbulence components can 
usually made over even with the rms meter. should pointed out 
that any tests operation this type instrument a-c fields, tanks, 
any situation where ground loops exist, will result high periodic noise 
levels such shown Fig. This noise level can usually eliminated 
drastically reduced changing the grounding configuration. 

example strip chart recording from turbulent tank flow shown 
Fig. Attempting sample every 0.01 sec yielded noisy result and the 
very poor fit reconstructed from the series had average point error about 
40% the amplitudes. Digitizing the short-term average recording like 
that shown Fig. yielded the Fourier components shown Table The 


= 


—_ 


FIG. 8.—Development Surface Turbulence IHNC Lock Outflow 


a = 
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fit was considered good showing average error 3%. The high amplitude 
components are seen below and almost energy contained 
above Hz. Fig. plot this frequency spectrum. 

can concluded from the laboratory tests considered, that this instrument 
capable detecting low frequency velocity spectral distributions natural 
turbulent water flow. Although the manual analysis methods used here are 
laborious and would not practical for large volumes data, tape recording 
and automated methods will yield spectral information for significant number 
data stations. Turbulence intensity can determined either the use 


TABLE 3.—Examples Measurements Made IHNC Lock with rms Voltmeters; 
Was Determined over 5-sec Time Average 


rms fluctuating 
Probe depth, Time average velocity, velocity, feet 
feet feet per second per second 
(meters) (meters per second) (meters per second) 


(1) (4) 


near surface 


Turbulence 
intensity, 


(0.61) 
(1.22) 


(1.83) 


FIG. 9.—Example Recording Longitudinal Component Turbulent Flow from 
Lock Outflow 


rms values from the digitized recordings rms voltmeter. previously 
covered, broad-band voltmeter measurements are used, noise considerations 
will important aspect data interpretation. 


Tests 


The field data collected included broad-band rms measurements well 
the strip chart recordings for the spectral analysis procedures The 


3.2 1.4 0.44 
(0.98) (0.43) 
2.9 0.3 0.10 
(0.88) (0.09) 
2.2 1.5 0.68 
(0.67) (0.45) 

1.8 1.2 0.67 
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sites will considered separately and, due the evolution the methods, 
the techniques and data are not the same each site. The differences are 
reviewed and the methods with results are compared the conclusion the 
end this section. 

IHNC.—The IHNC Lock the passage between the Mississippi River and 
the IHNC New Orleans. The outflow through duct and from head 
approx (4.57 m). Figs. and 8(b) are photographs spaced approx 
sec apart the control gates were opened showing the surface evolution 
the turbulent region. This flow can visually characterized large eddies 
with considerable entrained air the surface. The flow contained extensive 
apparently random mixing the center which was characteristic the radial 


TABLE Fourier Coefficients Longitudinal Turbulence Data Re- 
corded IHNC Lock Outflow 


Frequency, 
Hertz 


outward flow developed. Fronts turbulent water arrived the probe several 
seconds apart yielding very low frequency component addition those 
which appear the analysis covered. Considerable correlation between the 
magnitude the longitudinal and transverse velocities could observed within 
the flow they are indicated the instrumentation. 

Voltmeter measurements were made for several probe depths and the rms 
velocities were observed over approx sec. Examples the data are shown 
Table which taken the magnitude the vector sum the 
transverse channel and longitudinal channel velocities. The extensive 
randomness the flow made the designation specific direction transverse 
virtually meaningless. The values turbulence intensity the last column 


(1) (2) (3) 
5.97 0.00 0.0 
2.58 0.74 0.5 
1.0 
0.08 0.55 1.5 
1.73 —0.08 2.0 
0.88 2.5 
0.11 3.0 
0.27 3.5 
4.0 
0.52 4.5 
5.0 
5.5 
0.72 6.0 
6.5 
0.06 7.0 
0.28 —0.40 7.5 
0.16 0.38 8.0 
0.34 8.5 
0.18 9.0 
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indicate high levels turbulence compared maximum values about 
0.43 found McQuivey (4) the Columbia River Estuary. The data point 
depth was included show the approximate lower limit for the 
rms noise relative high flow. This measurement was made during quiet 


FREQUENCY 


FIG. 10.—Average Relative Amplitude Versus Frequency Longitudinal Turbulence 
IHNC Lock Outflow 


SEC 


FIG. 11.—Recording Longitudinal Component Turbulent Flow from Mississippi 
River New Orleans 


period between turbulent fronts. The flow this site general highly turbulent 
and may represent upper standard for comparison with other flows. 

Fig. shows example strip chart recording made this site, and 
Fourier analysis the short-term average yielded the components tabulated 
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TABLE 5.—Example Fourier Coefficients Longitudinal Turbulence Data Missis- 
sippi River 


Frequency, 
Hertz 


Note: Remaining terms are less than 0.50. 


RELATIVE AMPLITUDE 


FREQUENCY 


FIG. 12.—Average Relative Amplitude Versus Frequency Longitudinal Turbulence 
the Mississippi River 


HY3 321 
(1) (2) (3) 
0.86 0.33 
0.43 —0.28 0.66 
3.71 3.11 1.00 
0.48 1.33 
0.04 1.66 
0.57 0.74 2.0 
-0.11 2.33 
0.62 2.66 
0.48 0.32 3.00 
0.26 0.32 3.33 
0.61 3.66 
0.51 0.27 4.00 
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FIG. 13.—Recording Longitudinal Component Turbulent Fiow Mississippi 
River Approx 1,000 (304 Downstream from Working Tugboat 


FREQUENCY 


FIG. 14.—Relative Amplitude Versus Frequency Longitudinal Turbulence Recorded 
Mississippi River Downstream from Operating Tug 
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Table The most intense component 0.5 with more high amplitudes 
being found below 3.5 Hz. The mean error the fit this recording curves 
reconstructed from the Fourier components 3%. plot the average 
two these data increments shown Fig. indicating the concentration 
energy frequencies below Hz. 


TABLE 6.—Measurements and for Flow Under Seabrook Bridge 


rms 
Velocity, 
feet 
per second feet per per second 
(meters per (meters per 

Depth, second) (meters second) 

(meters) channel channel second) channel channel 


(2) (3) (5) (6) (7) 
(a) Moderate Light Turbulence 


Velocity, 
feet 


approx 
3.0 1.5 
(0.1) (0.46) (0.46) 


Mississippi River.—An example 3-sec data segment from recordings taken 
the Mississippi River seen Fig. 11. The probe was positioned (2.44 
deep and approx (3.05 from barges moored near Prytania Street 
New Orleans. Visual assessment the turbulent character the flow 
not usually possible the river because the wind ripple the surface. 
High noise levels were experienced due potentials induced the water around 


1.0 4.5 3.8 5.9 0.3° 0.06 
(0.30) (1.37) (1.16) (1.80) (0.12) (0.09) 
2.5 0.8 3.1 0.4 0.4 0.13 
(0.76) (0.91) (0.24) (0.94) (0.12) (0.12) 
4.0 0.0 2.5 1.5 1.5 0.60 
(1.22) (0.76) (0.46) (0.46) 
(5) Light Turbulent Flow 
1.0 2.4 4.5 5.1 0.25 0.20 0.04 
(0.30) (0.73) (1.37) (1.55) (0.08) (0.06) 
(c) Ship Passage 
4.0 0.5 3.1 3.0 3.0 0.97 
(1.22) (0.91) (0.15) (0.94) (0.91) (0.91) 
(d) Apparent Strong Turbulence 
(0.65) (0.21) (0.27) 
max 1.5 1.5 0.71 
(0.46) (0.46) 
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the barges which are wired with 60-Hz power. Table tabulation the 
Fourier components. Notice that the energy again concentrated the very 
low frequency region, with the higher amplitude components below Hz. Fig. 
plot relative amplitude versus frequency resulting from the average 
two spectra. The peak about 1.5 seemed characteristic data from 
the river, probably due boundary patterns generated the relatively periodic 
barge assembly which extended river from the observation site. The steady- 
state flow during these observations was 1.1 fps (0.33 m/s) and the rms measured 
the voltmeter was 0.20 rms fps (0.06 m/s) which yields value about 
0.18. McQuivey (4) measures values longitudinal component turbulence 
intensity far from the bank the Mississippi River Vicksburg high 
0.206 close the bottom and low 0.055 near the surface. The measurement 
made during this investigation was very close bounding surface, and compares 
well with McQuivey’s data from near the bottom (4). 


TABLE 7.—Summary Turbulence Field Data 


Visual Average 
indication turbulent 
intensity Spectral 
Location turbulence level characteristics 


(1) (2) (3) (4) 


IHNC High turbulence Chaotic with most the 
energy below 
Mississippi River Moderate turbulence Energy below with 
near barge peaks 1.0 Hz, 1.66 Hz, 
and 3.66 
Mississippi River Moderate turbulence Energy below with 
downstream high peak 1.0 
from operating 
tug 
Seabrook Bridge Light turbulent flow Spectrum not recorded 
Seabrook Bridge Light moderate Spectrum not recorded 
turbulence 
Seabrook Bridge High turbulence Spectrum not recorded 
Seabrook Bridge Ship passage Spectrum not recorded 


Fig. shows recording made the same location the Mississippi River 
except the wheel wash from tug boat operating about 1,000 (304 
was propagating across the probe. comparison this recording with Fig. 
shows considerable increase the relative amplitudes the lower frequency 
content the turbulence. plot the intensity function frequency 
shown Fig. 14. The relative concentration energy below Hzis significantly 
higher for generated the boat. The average error the fit 
the Fourier analyzed functions was less than 5%. 

Seabrook Bridge.—Flows under the Seabrook Bridge can representative 
highly turbulent and nonturbulent conditions depending wind, tide, and 
the location the probe. Voltmeter techniques were used this location during 
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strong flow period. These conditions could not found again for strip recording 
and spectral analysis. Table shows data and calculation turbulence intensity 
along with general visual description the flow for several probe positions 
indicated. 

The data for moderate turbulence were recorded flow between corner 
bounding wall and obstructions consisting old piles. The analysis 
function depth therefore not useful concept here, but several things 
can learned about the instrument from this data. The values turbulence 
intensity are within the reasonable limits for slightly turbulent flow and moderate 
turbulence. The region least turbulence yielded data labeled 
and the intensity considerably smaller than indicated the position with 
moderate turbulence. During the measurements barges under tow passed through 
the channel resulting increased turbulent character the flow. The value 
0.97 for the intensity quite high and even exceeds that found the IHNC 
Lock. This, however, transient condition, and these values represent the 
extreme. 

sequence measurements was made the flow where strong turbulence 
was visible the form mixing and eddies. The average lower values 
turbulence intensity were found about 0.4 and the higher values averaged 
near 0.7. These values seem bracket those obtained for the more turbulent 
flow regions measured McQuivey (4). 

Subsequent data were collected from the low drawbridge across the channel 
with flow ranging between 0.4 fps and 1.0 fps (0.12 m/s and 0.30 m/s). 
The rms measurements indicated between 0.1 fps rms and 0.2 fps rms which 
just above the expected noise levels flow. 


The summary turbulence intensity results tabulated Table shows that 
the level turbulence observed the different field test sites has been quantified 
and comparisons turbulent flows can made. The spectral plots the 
turbulence the river show definite periodic components while more chaotic 
spectra were observed the IHNC Lock. concluded that with this instrument 
method has been developed quantitatively compare turbulence intensity 
and obtain spectral signatures turbulent flows the field. 

Investigations the flows between Lake Ponchartrain and its connections 
the Gulf Mexico will conducted using these methods and results published 
later date. Automated data recording and digitizing will facilitate correlation 
calculations from which the characteristic mixing lengths can obtained. 
Scaling these data for comparison laboratory experiments such those 
done Etheridge and Kemp (2) will also attempted later date. 


Buchhave, P., George, K., Lumley, L., Turbulence With 
the Laser Doppler Annual Review Fluid Mechanics, Van Dyke, 
Wehausen, Lumley, eds., Vol. 11, pp. 443-503, 1979. 

Etheridge, W., and Kemp, H., Turbulent Flow Downstream 

Rearward Facing Journal Fluid Mechanics, Vol. 86, pp. 545-566. 
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McQuivey, S., and Measuring Techniques Turbulence Characteristics 
Open Channel USGS Professional Paper 802-A, United States Geological 
Survey, 1973. 

McQuivey, S., Turbulence Data From Rivers, Conveyance Channels, 


and Laboratory USGS Professional Paper 802-B, United States Geological 
Survey, 1973. 


The following symbols are used this paper: 


amplitude coefficients for cosine terms Fourier expansions; 
amplitude coefficients for sine Fourier expansions; 
time average resultant velocity 


average rms fluctuations for transverse and longitudinal velocity; 
output function for data set; 

index numbering terms Fourier expansion; 

time, seconds; 

extended time average longitudinal velocity; 


rms value fluctuating longitudinal velocity; 
extended time average transverse velocity; 


rms value fluctuating transverse velocity; 
instrument channel (x); and 
instrument channel 
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ULTIMATE DIMENSIONS LOCAL 


and George Members, ASCE 


The purpose this paper describe analytical method estimating 
the ultimate depth scour that based experimental observations made 
during the early stage the scouring process. Because scour approaches its 
limit asymptotically, current practice specify practical scour limit. This 
practical scour limit based personal judgment and can vary widely among 
individuals. The method propose defines scour limit independent individual 
judgment. 

Laursen has written (see Ref. pp. 180-181; see also Ref. 48): 


the following general characteristics which should basic any detailed analysis 
local scour [are]: 


The rate scour will equal the difference between the capacity for transport 
out the scoured area and the rate supply material that area. 

The rate scour will decrease the flow section enlarged. 

There will limiting extent scour. 

This limit will approached 


Among the persons concerned with local scour, there seems general 
agreement that Laursen’s characteristics describe the scour process. Herein 
the writers will take close look how scour approaches the limit with respect 
time, the aforementioned characteristic and suggest method analysis 


Hydr. Engr., United States Dept. Agr., Sci. and Education Administration, 
Agricultural Research, St. Anthony Falls Hydr. Lab., Third Avenue Southeast Mississippi 
River, Minneapolis, Minn. 55414. 

Engr., United States Dept. Agr., Sci. and Education Administration, Agricul- 
tural Research, St. Anthony Falls Hydr. Lab., Third Avenue Southeast Mississippi 
River, Minneapolis, Minn. 55414. 

Hydr. Engr., United States Dept. Agr., Sci. and Education Administration, 
Agricultural Research, St. Anthony Falls Hydr. Lab., Third Avenue Southeast Mississippi 
River, Minneapolis, Minn. 55414. 

Note.—Discussion open until August 1981. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication July 1980. 
This paper part the Journal the Hydraulics Division, Proceedings the American 
Society Civil Engineers, Vol. 107, No. HY3, March, 1981. 
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that will permit evaluation characteristic the rate scour, and characteristic 
the limit local scour. For this examination, assume the rate supply 
the scoured area zero, that the rate scour for characteristic 
will involve only the capacity for transport out the scoured area. 


Researchers have found that the shape local scour independent time 
suitable normalization parameters are selected nondimensionalize the units 
measurement (see Ref. pp. 183-184; Ref. 1197; Ref. pp. 283-287; 
and Ref. 50). However, the progression the scour with time has not 
been satisfactorily defined over the entire time range. This shown Anderson 
(Ref. 52, Fig. 2.17), who has plotted the scour depth against the logarithm 
time for the Laursen and Tarapore data. Regarding his Fig. 2.17, Anderson 
states (Ref. pp. 51-52): 


The curves show that the extent the scour increases linearly over the 
range shown with the logarithm the time. The linear relationship 
implies that the extent the scour increased indefinitely, contrary 
the hypothesis stated Laursen. must remembered, however, 
that such logarithmic function empirical and based fitting the 
relationship experimental points. fitting such curve, was often 
necessary omit points near the origin Fig. 2.17, because obvious 
that the logarithmic function cannot hold for small values time. 
the beginning the experiment, when [the time] zero, [the relative 
scour depth] also equal zero rather than minus infinity required 
logarithmic function. For larger values [the time], appears 
that the logarithmic function fits the data better, while for intermediate 
range, either logarithmic exponential function might suitable. 
Although experiments not indicate the condition clearly, 
plausible that for large values time, when the local shear stress 
approaching the critical value for the sediment, the scour-time relation 
may again depart from the logarithmic form. the absence analytical 
solution that describes the entire range, recourse must had empirical 
descriptions the various significant parts the entire curve. 

Note that, virtue the logarithmic character the development 
the scour region with time, practical equilibrium [emphasis supplied] 
reached after relatively short time, after which the increase the 
depth and extent [sic] scour becomes virtually imperceptible. 


The definition practical equilibrium can interpreted differently each 
person. When practical equilibrium reached? When does scour stop, either 
practically totally? agree that there are probably undefined practical 
limits the progression scour, but the random nature turbulent flow 
eddy bursts suggests that scour never stops completely. one experiment, 
the results which will presented later, the writers noted that scour continued 
even after 14-1/4 months. True, one might have look for several minutes 
detect the movement sediment particle, but movement was still continuing. 
The analytical method will describe determines the ultimate dimensions 
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scour—the asymptotic limit recognized most researchers—and thus removes 
the vague limit suggested the term practical equilibrium. 


Linear 


The writers’ ideas originated with and are extension the analytical method 
suggested Thomas (5,6,7) forecast the time development scour, and 
extrapolate experimentally-obtained values scour considerably longer 
times than are practical experiments. Thomas developed theoretical equation 
for the velocity scour—the scour depth divided the time since the beginning 
scour. His theoretical equation, presented Refs. (Eq. 75), and (Eq. 


Venax 
cosh arctanh 

Van 


cSpv, 


which the theoretical velocity which the scour depth increases; 
the theoretical scour depth; the time; v,, numerically, the fall 
velocity the sediment grain still fluid; the density the fluid; 
the density the sediment; the acceleration gravity; the 
velocity the jet the unscoured bed level; coefficient hydrodynamic 
aerodynamic resistance the sediment; the shaded area the sediment 
grain; and the volume the sediment grain. 

Thomas has computed theoretical values and has plotted log against 
log Fig. reproduced from Thomas’ paper (Ref. Fig. 5). Note that 


time 


FIG. 1.—Logarithmic Relationship Between Theoretical Velocity Scour and Time 
(Ref. Fig. 


"nth one 
2 
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for short time periods—less than 0.01 sec for the water data depending 
the sand size, and 0.5 sec for the air data—the velocity scour is, according 
Thomas, independent time. check Thomas’ statement that the velocity 
scour independent time, the writers differentiated the logarithmic form 
Eq. with respect the logarithm time for v,,, and evaluated 
the differential, various times. The evaluation shows that, 
although the slope negligible for small values the time, theoretically the 
velocity scour not independent time. (Experimental data presented 
Thomas show that this time independence may vary with scale, some data 
suggesting period long 50,000 sec—14 h.) (See Ref. Fig. for full-scale 
tests truss piers.) 

With increasing time and after transition, the logarithmic curve for Eq. 
is, also according Thomas, straight line again, but with average slope 
—0.948 for all the lines drawn Fig. Included are all sand sizes and 
both air and water. The equation this line has the form 


which and coefficients and exponents, respectively. 

noted previously, the average theoretical value for the curves 
Fig. given Thomas 0.948. However, the writers’ evaluation the 
differential the logarithmic form Eq. gave theoretical values 
1.000 for times exceeding 0.1 sec for 0.184-cm material. This value close 
agreement with the Thomas’ curves shown Fig. 

scour constant for all times scour. This constancy not supported 

Experimental data presented Thomas (Ref. Fig. for several structural 
forms that induce local scour have 0.931. Using this value Eq. 
becomes 


which the experimental scour depth. Thus, according the theory 
and experimental results presented Thomas, after the initial scour period 
the scour increases slowly but without limit. obtain practical equilibrium 
practical limit the scour for the Thomas method analysis, is, 
for the linear semilogarithmic method, also necessary specify the time period 
the scour progress. Therefore, this method analysis has the same limitation 
that the linear semilogarithmic method examined the preceding section. 


inspection the logarithmic plots experimental data presented 
Thomas—the plot Fig. (See Ref. Figs. and typical—shows that 
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the transition region much longer than suggested the theoretical plots. 
fact, the experimental data not exhibit the region where the scour velocity 
independent time. This not surprising: the region where independent 
time very short—this time period measured seconds fractions 
seconds. For these extremely short time periods, obtain meaningful 
experimental scour data very difficult impossible. The writers’ experience 
with the progression scour cantilevered pipe culvert outlets that 


satisfactory scour data cannot obtained for time periods less than about 
600 sec (10 min). 


Scour surface (a) 


Center rotation 


the plate 


FIG. Relationship Between Experimental Velocity Scour and Time: 
(a) Longitudinal Section; (b) Plot Experimental Data (Ref. Figs. and 


Many the experimental curves presented Thomas for the progression 
scour resulting from local disturbances caused baffles, piers, and spillway 
aprons exhibit the typical shape shown Fig. 

define the relationship between the logarithm the average velocity 
scour and the logarithm time, the writers looked for mathematical curve 
having shape similar the experimental curves Thomas’ plots log 
versus log and having asymptote from which the ultimate depth scour 
could computed. Such curve limb vertically-oriented hyperbola 
with the origin the hyperbola offset the y-axis. The general form 
the equation for the hyperbola shown Fig. 
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which the semitransverse axis; and the semiconjugate axis 
the hyperbola. have defined and dimensionless units when applying 
Eq. data from our research scour cantilevered pipe spillway culvert 
outlets. Our definitions and are: 


m 
log 


P P 


P 


log 


log 

P 
which Z,, the maximum depth scour measured from the tailwater surface; 
and the jet velocity and diameter the point where the jet plunges 
into the tailwater surface; and the time from the beginning scour. Because 


log V,t/D, appears both axes, the hyperbola rectangular, Eq. 
becomes 


HY3 
D, 
The writers have evaluated assuming values for and using experimental 
Asymptote 
-7 
FIG. 3.—General Form Hyperbola 
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values and compute The best value that value for which 
the standard error estimate, minimum. The procedure used is: 


Plot experimental values and and estimate the asymptote. Extrapolate 
the asymptote The distance the y-axis between the origin and 

Use the data compute values should constant.) 

Compute the standard error, 

Assume another value and repeat steps and 

Repeat step until the minimum value found. (The best fit values 
and are those for which minimum.) 


Fig. typical plot data obtained the writers. There will 
noted that the data, which covers scour periods ranging from 600 sec (log 


FIG. 4.—Typical Hyperbolic Logarithmic Plot Cantilevered Pipe Spillway Scour 


(2.31 cm)] 


4.60) 1.90 10’ sec (log 9.10) 0.5 increments the 
logarithm plus scour period 3.74 10’ sec (log V,t/D, 9.40), are far 


the right and approach the asymptote the hyperbola. The equation 
the hyperbola best fit 


and the ultimate asymptotic value the scour depth 


P 


(Experimental data for the 3.74 10” sec scour period did not become available 
until after Eq. had been evaluated. However, this omission affects only 
the numerical values Eqs. 12-15. The ommision the last data point does 
not affect the sense the hyperbolic logarithimic velocity-of-scour analytical 
method the comparison methods presented this paper.) 
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Linear Logarithmic inspection Fig. shows that the hyperbolic 
curve that best fits the data approximates straight line. The least squares 
fit straight line has the equation 


V,t (1—0.987643) Vit 0.012357 


P 
with correlation coefficient 0.999,970. This the form proposed Thomas 
(See Eq. 5). 

The time predicted Eq. reach the ultimate depth computed from 
the hyperbolic equation can computed substituting Z,,/D, computed from 
Eq. into Eq. 14. This time turns out V,t/D, or, for 
the data plotted Fig. which V,/D, 66.5 30,644 yr. 

Linear Semilogarithmic Method.—A similar comparison can made using 
the linear semilogarithmic method. The data plotted Fig. are replotted 
Fig. using semilogarithmic coordinates. Only after about 5-1/4 (log 


7 


FIG. 5.—Typical Linear Semilogarithmic Plot Cantilevered Pipe Spillway Scour 
(2.31 


the data fall approximate straight line. This agrees with Anderson’s 
Statement that the logarithmic function cannot hold for small values 
Therefore, only the data for scour periods greater than log (V,t/D,) 
were used develop the equation for the relationship between scour depth 
and time. The least squares fit equation 


5.5326 0.0486 log 


with correlation coefficient 0.850. 

The time predicted Eq. reach the ultimate depth computed the 
hyperbolic logarithmic equation V,t/D, 6.0 10” 2.9 10” yr. 

Summary Table.—The results obtained using the three computation methods 
are compared Table The slightly higher correlation coefficient for the 
hyperbolic logarithmic method indicates that the fit the data the linear 
logarithmic curve slightly inferior the fit the data the hyperbolic 
logarithmic curve. The correlation coefficient for the linear semilogarithmic curve 
relatively low. However, this low correlation coefficient for the linear 


6.5 
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semilogarithmic method and the high correlation coefficients for the logarithmic 
methods may mathematical and statistical rather than real. Another possible 
reason for the lower correlation coefficient for the linear logarithmic method 
that the actual scour depth was used for the semilogarithmic method, whereas 
the velocity scouring was used for the logarithmic methods. 


TABLE 1.—Comparison Methods 


Equation 


Hyperbolic Linear Linear 
logarithmic logarithmic semilogarithmic 
(2) (3) (4) 


(a) Comparison Equations 
Correlation 0.999,986 0.999,970 


Time reach ultimate 


scour depth --years 


(b) Time Comparison Scour Depths 


(8.76 hours) 
(3.65 days) 
(36.5 days) 


year) 
(10 years) 
(100 years) 

(1,000 years) 


(10,000 years) 
(100,000 years) 
(1,000,000 years) 
years) 


The ultimate scour depth, 6.98 computed the hyperbolic 
logarithmic method. 


Comparing the depths scour for the time periods shown Table 1(b) 
will permit each person determine what feels practical limit. 


The differences the predicted dimensions the scour between the various 
methods may small. The claim for the hyperbolic method that: (1) 
Determination the ultimate scour dimensions does not depend deciding 
what practical time; and (2) the limiting extent scour, Laursen’s third 
characteristic, can determined. the other hand, practical scour dimension 
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Item 
(1) 
Depth scour after 10” years 
5.69 5.64 5.84 
5.85 5.81 5.89 
5.97 5.97 5.94 
6.07 6.15 5.99 
6.16 6.32 6.03 
6.22 6.50 6.08 
6.28 6.69 6.13 
6.33 6.89 6.18 
6.38 7.08 6.23 
6.41 7.29 6.28 
6.98 bad 
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desired, the hyperbolic method analysis will provide prediction equation 
for any practical time period and also evaluate the difference between the practical 
and ultimate dimensions. 

The rate change the scour velocity, represented the slope the 
hyperbola relative its asymptote, maximum the beginning scour, 
decreases rapidly during the early stages scour, and becomes small the 
hyperbola approaches its asymptote. This change the scour rate with time 
agrees with Laursen’s second and fourth characteristics. 

The only scour dimension examined this paper has been the depth 
scour. However, view the comments made others (Ref. pp. 183-184; 
Ref. 1197; Ref. pp. 283-287; Ref. 50), and our unreported experiences, 
the comments made herein apply equally other dimensions local scour, 
the length and width the scoured area, and the location the scour 
with respect the pipe exit, pier, other source the local scour. 
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The following symbols are used this paper: 


semitransverse axis hyperbola; 

coefficient exponential equation; 

semiconjugate axis hyperbola; 

exponent exponential equation; height plate above original sand 
bed for Thomas’ air tests; 

coefficient hydrodynamic aerodynamic resistance sediment; 
pipe diameter; 

jet diameter point where jet plunges into tailwater surface; 
Froude number sediment grain for Thomas’ air tests; 
acceleration gravity; 

height approach channel for Thomas’ air tests; 

scour depth for Thomas’ air tests; 

discharge; 

shaded area sediment grain; 

time; 

volume sediment grain; 

jet velocity point where jet plunges into tailwater surface; 
numerically equal fall velocity sediment grain still fluid; 
velocity with which scour depth increases; 

velocity jet unscoured bed level; 

abscissa; 

ordinate; 

offset origin y-axis; 

maximum depth scour measured from tailwater surface; 
height pipe invert above tailwater level; 

angle deflector plate with horizontal for Thomas’ air tests; 
density fluid; 

density sediment; and 

standard error estimate. 


D, 
V, 


Subscripts 
experimental; and 
theoretical. 
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THREE-PARAMETER PROBABILITY DISTRIBUTIONS 
Donthamsetti Veerabhadra ASCE 


INTRODUCTION 


The United States Water Resources Council, while recommending the Pearson 
type distribution (P) with log transformation the data (log Pearson type 
distribution, LP) the base method for flood flow frequency analysis, identified 
continued study alternative distributions topic for future research 
(4). The principal consideration the use three-parameter probability distribu- 
tions like their flexibility that three sample statistics (usually 
the mean, variance, and skewness coefficient, are fitted the distribution when 
the method moments used for parameter estimation). comparative study 
four well-known three-parameter probability distributions the subject 
this paper. The distributions selected are: Lognormal (LN), Weibull (W), 
and LP. 
The three-parameter probability distributions are formed basically adding 
shift location parameter, their two-parameter counterparts. The location 
parameter (or function it, the case LP) serves bound (lower 
upper) for the random variable being fitted. Parameter not constant 
for given three-parameter probability distribution, but assumes different 
value for each combination the mean, variance, and skewness coefficient 
the random variable. Moreover, the bound becomes negative certain 
instances for LN, and This parameter has significant physical interpreta- 
tion when applied hydrologic variates such streamflow, precipitation, etc. 
These hydrologic variates cannot negative (unless made negative some 
transformation) one hand, and the other, imposition upper bound 
implies that there certain value (like flood which 
not physically meaningful. However, probability distributions with negative lower 
tails like the Gumbel were used for long time flood flow frequency analysis. 
such cases, since the quantiles the upper tail the distribution are 
interest, the fact that portion the distribution lies the negative region 
St. Johns River Water Management Dist., P.O. Box 1429, Palatka, Fla. 32077. 
Note.—Discussion open until August 1981. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications 
ASCE. Manuscript was submitted for review for possible publication February 26, 
1980. This paper part the Journal the Hydraulics Division, Proceedings the 
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the variate will not readily unless complete investigation 
the distribution made. Nevertheless, the should make 
sure least, that the area distribution the negative region (henceforth 
called area’’ for simplicity) the variate insignificant because 
the statistical parameters (SP) like are formed from the population 
variate values both positive and negative, whereas, the sample values are all 
positive. For the four three-parameter probability distributions selected, the 
bounds and the negative areas the distribution and selected quantiles are 
presented generalized fashion for comparison. Some guidelines are provided 
for selecting the best applicable three-parameter probability distribution for 
given hydrologic sample. 


For the four distributions selected, the probability density function, f(x), 
the relations the distribution parameters, and the SP, i.e., 


1.0 2.0 


and y,, and brief description, are follows. 


Lognormal: 


b? 


a) 2.0 
| b) Y= -2.0 


PROBABILITY DISTRIBUTIONS 


lognormal (LN2), with its origin shifted value, [see Fig. When 
negative, let represent the new location parameter; then, 
and (c’ distributed LN2 the reverse direction (mirror image 
the positive skewed distribution), with its origin shifted [see Fig. 
Eq. The unimodal and bell-shaped (with skew). Further 
details may found Aitchison and Brown (1). 

solve for the parameters and from given set p,, and 
the value determined first from Eq. trial. Then, and 
are solved explicitly. 

Weibull: 


The Weibull distribution, with modified form parameter equivalent 
the Gumbel’s Type III asymptotic extreme value distribution smallest 
values. Parameter the lower bound the distribution for both positive 
and negative y,. 

However, has lower limit —1.1396 (5). The Weibull distribution 
unimodal and bell-shaped when and has reverse J-shape 
given SP, parameter solved first, trial, from Eq. 2d, and then and 
are solved explicitly. 

Pearson Type (three-parameter gamma) 
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gamma (GA2) with its origin shifted LN, when negative, 
distributed GA2 the reverse direction. The and become 
equivalent when Like also unimodal and bell-shaped when 

Log Pearson Type 


exp (ac) 


a a a a 


x 


either positive negative depending the values and and 
exp (c). The many forms taken are analyzed length Bobee 
(2). Solution parameters from given available Ref. 


the method moments, parameters distributions are evaluated equating 
population SP, i.e., and the corresponding sample estimates. Thus, 
will more comprehensive comparison different properties 
three-parameter probability distributions made the basis SP. the 
following sections some properties the four distributions presented are 
evaluated terms the dimensionless variate Such transformation 
the random variable, meaningfully scales data, giving value unity 
for the mean i.e., 1.0. The transformed data, X,/X 
(in which the ith observed value and sample mean), being ratios 
the mean, afford easy for the results obtained. The standard deviation 
1.0, the distributions yield different results terms All results may 
classified and compared the basis two SP, and y,. The value 
the variance and represents the squared value coefficient 
variation data. For floods real space, Landwehr al. (6) evaluated 
the values mean coefficient variation (CV) for the hydrologic regions 
the United States for record lengths, yr, yr, and yr. Values 


a 

@ 
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varied from 0.39-0.99 for the regions and the standard deviations 
for are the range 0.15-0.43 (the extreme values being applicable 
the respective extreme values CV). These statistics indicate the general range 
encountered hydrologic samples. 

Ranges for and —9.0-9.0 for were chosen for investigation. 
However, the cases and LP, values less than were not investigated 
because cannot take —1.1396 and assumes U-shape for most 
negative values (7). For several combinations and chosen systemati- 
cally, the distribution parameters, and whence different results interest, 
were evaluated. 


The bound values different distributions are presented Tables 1-4. The 
and are bounded the lower tail when positive (Table 1), and 
they are bounded the upper tail when negative (Table 2). The 


TABLE 1.—Lower Bound Values (Expressed Ratios Mean) Lognormal and 
Pearson Type Distributions 


VARIANCE SKEWNESS COEFFICIENT ( 


0.50 1.00 1.50 2.00 3.00 4.00 5.00 7.00 9.00 
(2) (3) (4) (S) (6) (7) (8) (9) (10) 


LOGNORMAL DISTRIBUTION 


0.395 0.690 0.786 0.832 0.878 0.900 0.913 0.929 0.938 
-0.354/| 0.306 0.520 0.625 0.727 0.776 0.806 0.841 0.861 
-0.915| 0.018 0.322 0.469 0.613 0.684 0.726 0.775 0.804 
-1.708 -0.388 | 0.041 0.250 0.453 0.553 0.613 0.682 0.722 
-2.316 -0.700 -0.175 | 0.081 0.330 0.452 0.525 0.611 0.660 
-2.829 -0.963 -0.357 -0.061 | 0.227 0.368 0.452 0.550 0.607 
-3.281 -1.195 -0.517 -0.186 | 0.135 0.293 0.387 0.497 0.561 
-3.690 -1.404 -0.661 -0.300 [0.053 0.225 0.329 0.449 0.519 
-4.066 -1.597 -0.795 -0.404 -0.023 | 0.163 0.275 0.405 0.480 
~4.415 -1.776 -0.918 -0.501 -0.094 | 0.106 0.225 0.364 0.444 
-4.744 -1.945 -1.035 -0.592 -0.160 | 0.051 0.178 0.325 0.411 
-5.055 -2.104 -1.145 -0.678 -0.22 0.0004.0.134 0.289 0.379 
-7.562 -3.389 -2.033 -1.373 -0.729 -0.414 -0.225 -0.006 Lo.121 
-9.487 -4.376 -2.715 -1.906 -1.118 -0.732 -0.501 -0.232 -0.076 
11.109 -5.208 -3.290 -2.355 -1.446 -1.000 -0.733 -0.422 -0.242 
12.538 -5.940 -3.796 -2.751 -1.734 -1.236 -0.937 -0.590 -0.389 


PEARSON TYPE 3 DISTRIBUTION 


0.600 0.800 0.867 0.900 0.933 0.950 0.960 0.971 0.978 
0.106 0.553 0.702 0.776 0.851 0.888 0.911 0.936 0.950 
-0.265 |0.368 0.578 0.684 0.789 0.842 0.874 0.910 0.930 
-0.789 [0.106 0.404 0.553 0.702 0.776 0.821 0.872 0.901 
-1.191 -0.095 | 0.270 0.452 0.635 0.726 0.781 0.844 0.878 
-1.530 -0.265 | 0.157 0.368 0.578 0.684 0.747 0.819 0.859 
-1.828 -0.414 | 0.057 0.293 0.529 0.646 0.717 0.798 0.843 
-2.098 -0.549 0.225 0.484 0.613 0.6970 0.779 0.828 
-2.347 -0.673 -0.116] 0.163 0.442 0.582 0.665 0.761 0.814 
-2.578 -0.789 -0.193 | 0.106 0.404 0.553 0.642 0.744 0.801 
-2.795 -0.897 -0.265 | 0.051 0.368 0.526 0.621 0.729 0.789 
-3.000 -1.000 -0.333 + 0.0004 0.333 0.500 0.600 0.714 0.778 
-4.657 -1.828 -0.886 -0.414 | 0.057, 0.293 0.434 0.596 0.686 
-5.928 -2.464 -1.309 -0.732 -0.155| 0.134 0.307 0.505 0.615 


-7.000 -3.000 -1.667 -1.000 -0.333 40.0004 0.200 0.429 0.556 
-7.944 -3.472 -1.981 -1.236 -0.491 -0.118! 0.106 0.361 0.503 
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always bounded the lower tail (Table 3). The the only distribution 
the four with lower bound The has upper bound when 
negative, and for some combinations and when positive 
(Table 4). The combinations and for which has lower and upper 
bounds are denoted Class and Class respectively, Table When 
the has upper bound, develops U-form increased given 
this occurs because the constrained (from spreading) from both ends 


TABLE 2.—Upper Bound Values (Expressed Ratios Mean) Negatively Skewed 
Lognormal and Pearson Type Distributions 


VARIANCE SKEWNESS COEFFICIENT ) 
2 
“0.50 -1.00 -1.50 -2.00 -3.00 -4.00 -5.00 -7.00 -9.00 


(2) (3) (4) (3S) (6) 


(7) (8) 


(9) (10) 


LOGNORMAL DISTRIBUTION 


1.214 1.168 1.122 1.100 
0.05 2-354 1.694 1.480 1.375 1.273 1.224 1.194 1.159 1.139 
0.10 2.915 1.982 1.678 1.531 1.387 1.316 1.274 1.225 1.196 
0.20 3.708 2.388 1.959 1.750 1.547 1.447 1.387 1.318 1.278 
0.30 4.316 2.700 2.175 1.919 1.670 1.548 1.475 1.389 1.340 
0.40 4.829 2.963 2.357 2.061 1.773 1.632 1.548 1.450 1.393 
0.50 5-281 3.195 2.517 2.186 1.865 1.707 1.613 1.503 1.439 
0.60 5.690 3.404 2.661 2.300 1.947 1.775 1.671 1.551 1.481 
0.70 6.066 3.597 2.795 2.404 2.023 1.837 1.725 1.595 1.520 
0.80 6-415 3.776 2.918 2.501 2.094 1.894 1.775 1.636 1.556 
0.90 6-744 3.945 3.035 2.592 2.160 1.949 1.822 1.675 1.589 
1.00 7-055 4.104 3.145 2.678 2.223 2.000 1.866 1.711 1.621 
2.00 9.562 5.389 4.033 3.373 2.729 2.414 2.225 2.006 1.879 
3.00 11.487 6.376 4.715 3.906 3.118 2.732 2.501 2.232 2.076 
4.00 13.109 7.208 5.290 4.355 3.446 3.000 2.733 2.422 2.242 
5-796 4.751 3.734 3.236 


PEARSON TYPE 3 DISTRIBUTION 


1.400 1.4200 1.133 1.100 1.067 1.050 1.040 1.029 1.022 
0.05 1.894 1.447 1.298 1.224 1.149 1.112 1.089 1.064 1.050 
0.10 2-265 1.632 1.422 1.316 1.211 1.158 1.126 1.090 1.070 
0.20 2.789 1.894 1.596 1.447 1.298 1.224 1.179 1.128 1.099 
0.30 3.191 2.095 1.730 1.548 1.365 1.274 1.219 1.156 1.122 
0.40 3.530 2-265 1.843 1.632 1.422 1.316 1.253 1.161 1.141 
0.50 3.828 2.414 1.943 1.707 1.471 1.354 1.283 1.202 1.157 
0.60 4.098 2.549 2.033 1.775 1.516 1.387 1.310 1.221 1.172 
0.70 4.347 2.673 2.116 1.837 1.558 1.418 1.335 1.239 1.186 
0.80 4.578 2.789 2.193 1.894 1.596 1.447 1.358 1.256 1.199 
0.90 4-795 2.897 2.265 1.949 1.632 1.474 1.379 1.271 1.211 
1.00 5-000 3.000 2.333 2.000 1.667 1.500 1.400 1.286 1.222 
2.00 6-657 3.828 2.886 2.414 1.943 1.707 1.566 1.404 1.314 
3.00 7-928 4.464 3.309 2.732 2.155 1.866 1.693 1.495 1.385 
4.00 7.000 5.000 3.667 3.000 2.333 2.000 1.800 1.571 1.444 
5.00 9-944 5-472 3.981 3.236 2.491 2.118 1.894 1.639 1.497 


whereas the dispersion data (variance) about the mean increased. These 
combinations and are indicated the letter Table Table 
also shows, footnotes, the combinations and which the 
assumes reverse J-shape; all other instances (of Table 4), the 

For hydrologic data, the bounds which may considered desirable are: 
lower bound and upper bound equal very large value like, 
100 (which represents flood value one hundred times the annual mean, 
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TABLE 3.—Lower Bound Values (Expressed Ratios Mean) and Negative Areas 
Weibull Distribution 


VARIANCE SKEWNESS COEFFICIENT ( %? 
2 
7% -1.00 -0.50 0.00 0.50 1.00 1.50 2.00 3.00 4.00 5.00 7.00 9.00 
(2) (3) (4) (6) (7) (8) (10) (11) (12) (13) 


LOWER BOUND VALUES 


-3.484 0.676 0.847 0.879 0.900 

-9.026 -0.413 | 0.275 0.658 0.730 0.776 
~13.179 -0.998 0.516 0.619 0.684 
-19.052 -1.825 0.316 0.461 0.553 
~23.559 -2.460 0.162 0.340 0.452 
-27.358 -2.996 0.032 0.238 0.368 
-~30.705 -3.467 -0.0 0.148 0.293 
~33.732 -3.894 -0.185 | 0.066 0.225 
-36.514 -4.286 -0.280 0.163 
~39.104 -4.651 -0.369 -0.078| 0.106 
~41.537 -4.993 -0.452 -0.144/ 0.051 
~43.838 -5.318 -0.530 -0.205 0.000 
~62.411 -7.934 -1.966 -1.164 -0.705 -0.414 +074 
-76+662 -9.942 ~2.633 -1.651 -1.088 -0.732 -0.315 
~88.676-11.635 -3.195 -2.061 -1.411 -1.000 -0.519 
~99 «261-13.126 -3.690 -2.422 -1.695 -1.236 -0.698 


‘NEGATIVE AREAS’» PERCENT 


0.000 

0.147 0.006 

0.895 0.349 0.000 

3.080 2.360 1.013 

5.246 4.765 3.425 0.821 

74163 6.996 5.997 3.363 

8.833 8.965 8.360 6.166 1.490 

10.292 10.689 10.453 8.781 4.543 

11.575 12.203 12.295 11.117 7.567 0.286 

12.715 13.542 13.921 13.182 10.310 3.785 

13.735 14.736 15.363 15.009 12.745 7.227 

14.655 15.807 16.650 16.631 14.897 10.292 
20.642 22.643 24.676 26.453 27.445 27.216 25.390 13.428 
23.898 26.252 28.764 31.230 33.192 34.324 34.469 31.330 20.527 
26.031 28.575 31.339 34.159 36.600 38.364 39.347 38.944 35.050 24.742 
27.572 30.232 33.150 36.183 38.904 41.029 42.466 43.382 41.883 37.530 


TABLE 4.—Upper and Lower Bound Values Log Pearson Type Distribution 


VARIANCE SKEWNESS COEFFICIENT ( y, ) 
“1.00 -0.50 0.00 0.50 1.00 1.50 2.00 3.00 4.00 5.00 7.00 9.00 
«> (2) (s) (6) (8) (9) CBD 


LOWER BOUND VALUES(CLASS A) 
0.01 1.158 1.277 1.951 0.739 0.835 0.878 0.918 0.937 0.948 
0.05 1.275? 1.442 1,958 14.260] 0.202 0.520 0.653 0.767 0.818 0.846 
0.10 1.328 1.512. 1.967 4.429 0.207 0.425 0.4616 0.697 0.742 
0.20 1.573? 1.983 3.208 15.875] #8 0.084 0.355 0.487 0.559 
0.30 1.996 2.944 7.030 0.153 0.401 
0.40 2.835 5.408 27.885 0.035 0.268 
0.50 2.779 4.761 14.172 0.001 0.163 
0.60 2.748! 4,423 10,294 69.920 0.086 
0.70 4.220 8.559 33.070 * 0.037 
0.80 4.088! 7,600 21.972 * 0.012 
0.90 3.996! 7.000 17.005 * 0.002 
1.00 3.930! 6.594 14.278 * 

3.00 u 17.1291 49.132 * 
4.00 U 14. 706! 33.122! 86.689 
5.00 u U 13.669! 27.498! 59.603! 

UPPER BOUND VALUES(CLASS B) 


NOTES: - 
FOR CLASS A THE UPPER BOUND VALUE = = 1 - THE DISTRIBUTION IS REVERSE J-SHAPED 
FOR CLASS B THE LOWER BOUND VALUE = 0 2 - THE DISTRIBUTION IS J-SHAPED 


= K > 100 U - THE DISTRIBUTION IS U-SHAPED 
0 < K < .0005 


ececcccccecccce 
eeccccccecccc 
ececccccecce 
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0.01 0.938 0.946 0.957 0.963 
0.05 0.860 0.880 0.904 0.918 
0.10 0.803 0.830 0.864 0.884 
0.20 0.721 0.760 0.807 0.836 
0.30 0.658 0.706 0.764 0.799 
0.40 0.605 0.660 0.727 0.767 
0.50 0.559 0.620 0.695 0.740 
0.60 0.516 0.584 0.466 0.715 7 
0.70 0.478 0.550 0.639 0.692 ) 
0.80 0.442 0.519 0.614 0.671 
0.90 0.408 0.490 0.591 0.651 
1.00 0.376 0.462 0.569 0.632 
2.00 0.117 0.240 0.390 0.480 
3.00 0.069 0.253 0.363 
4.00 -0.249 -0.075| 0.137 0.265 
5.00 -0.396 -0.202! 0.035 0.178 
0.01 
0.05 
0.10 
0.20 
0.30 
0.40 
0.50 
0.60 
0.70 
0.80 
0.90 
1.00 
2.00 
3.00 
4.00 
5.00 
0.960 0.9671 
0.877 0.893! 
0.789 0.814 
0.635 0.673 
0.499 0.550 
0.381 0.441 
0.281 0.347 
0.199 0.267 
0.133 0.200 
0.083 0.145 
0.047 0.101 
0.068 


desirable bounds lie along narrow band and values (see the dividing 


346 MARCH 1981 HY3 
which may never reached practice). Tables and show that these 


lines Tables and for and when positive, and for 
and LP, general. fact, the cases which the lower bound and the 
upper bound are equivalent their two-parameter counterparts for LN, 
and and LN2 for (note, however, that LN2 not the two-parameter 
probability distribution but only special case LP). two-parameter probability 


TABLE and Relations for Two-Parameter Probability Distributions 


(1) (2) (3) (4) 


0.0100 0.2000 
0.0500 


0.1000 0.9803 0.0245 
0.2000 1.4311 0.8944 0.4148 
0.3000 1.8075 1.0954 0.7020 
0.4000 1.2649 0.9410 
0.5000 1.1514 
0.6000 1.5492 1.3428 
0.7000 3.0956 1.6733 1.5204 
0.8000 3.3988 1.7889 1.6879 
1.8472 
3.0000 10.3923 3.4641 
4.0000 14.0000 4.0000 
17.8885 4.4721 6.6188 


distributions, and are related follows: 


LN2: +30, 


W2: Parameter and related as: 


GA2: Parameter and are related Eq. 2d: 


The numerical values and for two-parameter probability distributions 
are presented Table may seen from Tables and the 
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bounds the three-parameter probability distributions become less desirable 
value deviates more and more from the two-parameter probability distribu- 
tion value given value deviates more and more from the 
two-parameter probability distribution value given y,. This suggests that, 
least from the bounds point view, comparison sample skewness 
coefficient (CS) with two-parameter probability distributions the sample 
variance will indicate which three-parameter probability distribution more 
suitable for application. example, assume that and were calculated 
0.5 and 1.5, respectively, for sample. Since close GA2, 
will have better acceptable bounds than the other three distributions. 

When negative, the upper bounds LN, and appear generally 
low (Tables and 4); their suitability should determined conjunction with 
the type and the largest item data and the negative area the distribution. 


Areas 


The LN, and enter the region negative values for several combinations 
and y,. such cases, their application analyze hydrologic variables 


TABLE 6.—Negative Areas for Positively Skewed Lognormal and Pearson Type 
Distributions, percentage 


VARIANCE SKEWNESS COEFFICIENT Cy) 


2 
k 1.50 2.00 3.00 4.00 5.00 7.00 9.00 
(1) (4) (S) (6) (7) (8) (9) (10) 


LOGNORMAL DISTRIBUTION 


0.000 

0.001 

0.319 0.005 

1.762 0.384 0.002 

3.896 1.747 0.263 0.000 

6.184 3.770 1.368 0.103 

8.383 6.009 3.180 0.852 

10.412 8.223 5.312 2.359 0.000 

12.254 10.311 7.506 4.313 0.105 

13.920 12.238 9.630 6.435 0.798 

15.426 14.000 11.628 8.560 2.163 
24.984 25.285 24.821 23.672 19.804 14.488 8.376 0.000 
29.846 30.932 31.373 31.231 29.631 26.824 23.185 14.067 3.766 
32.875 34.391 35.319 35.715 35.313 33.860 31.735 26.096 18.949 
34.986 36.770 38.000 38.722 39.035 38.372 37.113 33.479 28.771 


PEARSON TYPE 3 DISTRIBUTION 


0.000 

0.260 

1.672 0.047 

3.837 1,060 

6.170 3.136 

8.412 5.578 0.354 
10.4746 8.013 2.670 
12.347 10.298 5.513 
14.035 12.390 8.303 
15.558 14.288 10.882 
25.177 26.091 26.414 25.390 
30.041 31.816 33.386 34.469 31.644 
33.064 35.277 37.424 39.347 41.113 
35.167 37.641 40.111 42.466 45.978 44.235 


0.05 
0.10 
0.20 
0.30 
0.40 
0.50 
0.60 
0.70 
0.80 
0.90 
1.00 
2.00 
3.00 
4.00 
5.00 
0.10 
0.20 
0.30 
0.40 
0.50 
0.60 
0.70 
0.80 
0.90 
1.00 
2.00 
3.00 
4.00 
5.00 
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(which are positive) may not meaningful. However, they may still applied 
when the variate the upper tail distribution (like flood flows) are 
interest the area negative lower tail insignificant. For the three 
distributions mentioned, the negative areas are presented percentages total 

The negative areas distributions may obtained evaluating the cumulative 
distribution function, between the limits and the lower (negative) 


TABLE 7.—Negative Areas for Negatively Skewed Lognormal and Pearson Type 
Distributions, percentage 


VAR TANCE SKEWNESS COEFFICIENT ( 
2 
-0.50 -1.00 -1.50 -2.00 -4.00 -5.00 -7.00 -9.00 
(1) (2) (3) (4) (S) (6) (7) (8) (9) (10) 


LOGNORMAL DISTRIBUTION 


0.000 0.000 0.002 0.007 0.027 0.049 0.067 0.091 
0.028 0.137 0.287 0.425 0.608 0.699 0.742 0.768 
0.410 0.837 1.177 1.400 1.608 1.662 1.658 1.598 
2.286 2.981 3.347 3.499 3.508 3.386 3.242 2.982 
4.565 5.164 5.366 5.358 5.109 4.806 4.531 4.094 
6.728 7.115 7.119 6.945 6.457 5.995 5.607 5.021 
8.666 8.822 8.634 8.310 7.611 7.012 6.528 5.816 
10.380 10.315 9.954 9%.497 8.615 7.898 7.332 6.511 
11.896 11.629 11.115 10.542 9.500 8.680 8.043 7.128 
13.244 12.797 12.146 11.471 10.290 9.380 8.681 7.682 
14.449 13.841 13.070 12.304 11.000 10.012 9.257 8.186 
15.534 14.782 13.904 13.058 11.645 10.586 97.782 8.646 
22.503 20.891 19.374 18.048 15.976 14.488 13.381 11.837 
26.199 24.197 22.386 20.837 18.454 16.760 15.505 13.757 
28.580 26.356 24.377 22.700 20.135 18.319 16.974 15.103 
30.278 27.912 25.824 24.063 21.379 19.482 18.079 16.124 
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0.000 0.006 0.000 0.002 0.013 0.037 0.069 
0,023 0.110 0.253 0.420 0.727 0.937 1.055 
0.392 0.814 1.221 1.557 1.984 2.151 2.162 
2-285 3.069 3.606 3.932 4.137 4.005 3.733 
4.598 5.352 5.764 5.926 5.787 S.350 4,832 
6.789 7.369 7.592 7.569 7.097 6.395 5.673 
8-747 9.116 9.144 8.944 8.175 7.244 6.352 
10.475 10.633 10.477 10.117 9.084 7.957 6.919 
12.000 11.961 11.638 11.133 9.869 8.570 7.405 
13.354 13.135 12.661 12.027 10.556 9.106 7.830 
14.564 14.181 13.571 12.821 11.166 9.581 8.206 
15.651 15.120 14.387 13.534 11.713 10.006 8.543 
22.616 21.161 19.654 18.139 15.260 12.773 10.737 
26.297 24.398 22.505 20.652 17.218 14.313 11.966 
28.665 26.503 24.376 22.313 18.527 15.351 12.799 
30.353 28.016 25.730 23.522 19.489 16.119 13.419 


bound. Only has equation for F(x) its closed form, given 


F(x) 


the cases and numerical integration methods were used. Note 
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0.01 0.105 
0.05 0.763 
0.10 1.526 
0.20 2.775 
0.30 3.772 
0.40 4.603 
0.50 5.316 
0.60 5.941 
0.70 6.498 
0.80 6.999 
0.90 7.455 
1.00 7.872 
2.00 10.795 
3.00 12.578 
4.00 13.839 
5.00 
0.01 0.134 0.181 
0.05 
0.10 1.985 1.740 
0.20 3.114 2.574 
0.30 3.870 
0.40 
0.50 4.890 
0.60 5.268 4,107 
0.70 5.590 
0.80 4.531 
0.90 
1.00 6.341 4,861 
2.00 7.793 5.882 
3.00 8.614 6.462 
4.00 9.175 6.861 
5.00 9.596 7.162 
@ 
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Tables and show that given negative areas the three distributions 
increase with indicating thereby that these distributions spread more and 
more into the negative region increased. When positive they 
enter the negative region when greater than the two-parameter probability 
distribution value. 


One the main uses probability distributions hydrology predict 
the magnitude some hydrologic variate associated with some exceedance /non- 
exceedance probability. These probabilities are more commonly expressed 
recurrence intervals return periods, when the data analyzed are annual 
events. For the value variate with nonexceedance probability equal 
given 1/F when data the lower tail distribution are concern, 
i.e., the annual smallest events (ASE) like low flows, and when 
data the upper tail distribution are concern, i.e., the annual largest 
events (ALE) like flood flows. 

The quantiles may evaluated from the tables normal distribution 
mean and standard deviation The quantiles may evaluated 
Eq. The tables available Ref. may used evaluate quantiles 
interest for and LP. These tables provide values the equation 


+o, K’ 


which Pearson quantile nonexceedance probability and and 
the mean and standard deviation respectively. The value 
and being common for both and LP. After determining parameters 
use Eq. the quantile, then, given exp (In Ref. the 
values this paper were called values which may noted the reader 
avoid confusion.) 

The values selected frequencies were evaluated for the four distributions 
and presented Tables and respectively, for the ASE and the ALE. 
Some the ASE are negative for LN, and because these distributions 
enter the negative range The reciprocal the negative area (Tables 
and gives for each these distributions the limiting (of ASE); when 
this exceeded, the quantiles will negative. 

Tables and show that for some combinations and y,, the values 
for given are practically equal for some all the four distributions. 
This indicates that under certain circumstances (shown quantitatively for different 
ranges and Tables and will not really matter which probability 
distribution chosen the method moments used for estimating parameters. 
Further, observed that the values different distributions for 
differ greatly when these distributions have properties which may not 
acceptable hydrologic analysis. example, take the case sample 
for which 1.0 and 1.0 and the data are flood flows. Table shows 
that, for this case, the value differs greatly from the other three 
distributions for 100 yr; the 500-yr value less 29%-35%. The 
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bounded upwards, and the bound value 3.93 (Table 4). The LN, 
and enter the negative range and the lower bounds are —2.104, 
—0.530, and —1.00, respectively, and the negative areas are 14%, 14.9%, and 
14.3%, respectively. These properties, general, may unacceptable for 
analyzing hydrologic samples. the case ALE, quantiles are, general, 

The results presented Tables 1-9 were based theoretical evaluation 
distributions. Some combinations and presented these tables, 
like low and high vice versa, may rarely never found 
practical examples hydrology. Fortunately, for these combinations 
and that the distribtuions have properties not acceptable hydrologic 
frequency. 


Practicat 


Different results presented the foregoing sections show that the skewness 
coefficient the two-parameter counterpart three-parameter probability 
distribution (LN2 the case LP) can serve index indicate the 
behavior the three-parameter probability distribution given SP. This 
will denoted y,, henceforth. The value y,, evaluated sample 
variance, When positive: (1) The and will enter the negative 
and (3) the lower bounds LN, and are greater than zero, 
bound greater than zero for both positive and negative CS. 
When negative, the LN, and are bounded upwards. Further, the 
more the deviates from y,,, the more undesirable the properties the 
three-parameter probability distributions become for hydrologic analysis. Thus, 
general, for given hydrologic sample, the three-parameter probability 
distributions, which has its y,, closer CS, will have more appealing 
properties than others. However, mentioned earlier, for certain values 
SP, some the quantiles (for given different distributions may not 
differ greatly and the choice distribution only academic interest. 
Applications these concepts are illustrated through the following examples. 

Annual flood flow data from two gaging stations (St. Johns River near DeLand, 
Fla., and Cedar River Cedar Rapids, were selected for analysis. The 
data were first converted into dimensionless variates, K,, the formula 
X,/X. The unbiased estimates sample variance and skewness coefficients 
and CS,) are calculated usually the following formulas: 


i=l 


which 1.0; and the sample size. 
However, Monte Carlo experiments conducted Wallis, al. (9) indicated 
that Eq. does not sufficiently correct for the sample bias the cases 
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LN, and Based the results Ref. Bobee and Robitaille (3) proposed 
the following formulas for the unbiased estimates skewness for LN, 
and distributions: 


Distribution: 


Distribution: 


Distribution: 


the aforementioned, mean the distribution the skewness coefficient 
for sample size Normally only one sample size available and 
substituted the biased estimate skewness coefficient given 


For the also, Eq. inadequate correct for the sample bias (8). 
The writer has evolved the method mixed moments for (8), 
which the skewness coefficient determined the basis 
sample logarithmic mean, instead using The value Yis calculated 


Since unbiased estimate, may regarded unbiased less 
biased than CS,. The LN, and are fit the two samples (by the 
method moments), first using the statistics given Eqs. 10-11, and 
then replacing CS, the new unbiased estimates, i.e., the case 
LN, CS,, the case etc. Table presents the estimates parameters 
and and the values for return periods, yr, 100 yr, 200 
yr, and 500 for the four distributions. For each distribution, the results obtained 
using the new unbiased estimate skewness are shown parentheses. 
analysis the results follows. 

Sample No. St. Johns River near DeLand, data are follows: 
Sample size (1934-1978); 8,126 cfs (230.13 0.18133; 
CS, 0.9041; 1.0745; 1.0083; CS, 1.0330; and CS,, 
1.3591. Skewness coefficients y,,, LN2, W2, and GA2 (at variance 0.18133) 
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are: 0.3526 and 0.8517, respectively (From 5-7, substituting 

Case 1.—Results obtained using CS, for all distributions: this case, 
enter the negative range will bounded upwards, and and 
will have lower bound greater than zero. However, since close 
the case will have better acceptable bounds. The has 
upper bound, exp (c) 11.450. The LN, and have lower bounds, 
—0.453, 0.313, and 0.058, respectively. The negative area 0.0045%, 
which very negligible. For yr, the quantiles and 
are practically equal, and values differ from and 0.3%-1.7% (Table 
10). Thus, although has better attributes, also gives equivalent results. 


TABLE 10.—Results Examples 


Values for Return Period 
a 


Distribution 


(1) 


Lognormal 
Weibull 
Pearson 


Log Pearson 


Lognormal 
Weibull 
Pearson 


Log Pearson 


Case obtained using the new unbiased estimates skewness 
for each distribution: this case, but CS,, and CS,, are 
greater than y,,. Therefore, will enter the negative range and 
and will have lower bound greater than zero. However, CS, practically 
equal and thus clearly the choice this case. With use 
the new unbiased estimates skewness, the statistical attributes 
and improved, but those and became less attractive. The LN, 
and have lower bounds, —0.236, 0.351 and 0.176, respectively, and 
has lower bound, exp (c) The negative area 0.0001%. 

Sample No. Cedar River Cedar Rapids, lowa.—The data are follows: 
Sample size (1903-1972). The writer had opportunity work this 


(a) St. Johns River near Deland, Fla. 

0.333 0.287 —0.453 2.061 2.266 2.468 2.732 
(0.156) (0.335) (—0.236) (2.089) (2.311) (2.533) (2.828) 

1.658 0.769 0.313 2.063 2.244 2.415 2.627 

(1.556) (0.722) (0.351) (2.085) (2.277) (2.458) (2.685) 

5.195 4.893 0.058 2.065 2.260 2.450 2.692 

(4.547) (3.749) (0.176) (2.088) (2.296) (2.498) (2.757) 

(1620.) (4.37 10°) (—269.7) (2.128) (2.379) (2.634) (2.981) 

(b) Cedar River Cedar Rapids, lowa 

0.837 0.262 2.563 2.855 3.142 3.515 

(0.662) (0.303) (—1.030) (2.581) (2.892) (3.200) (3.605) 

1.753 1.214 —0.081 2.561 2.819 3.060 3.359 

(1.639) (1.123) (—0.005) (2.577) (2.847) (3.102) (3.419) 

3.844 5.982 2.566 2.847 3.118 3.463 
(3.397) (4.565) (—0.344) (2.581) (2.875) (3.159) (3.524) 

4.963 1.460 2.570 2.798 2.795 3.217 
(9.988) (1.997) (2.655) (2.958) (3.241) (3.588) 
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sample different context the past, but did not get chance update 
data for the present study. 27,558 cfs (771.62 0.3957; 
0.8595; 0.9595; 0.9227; CS, 0.9360; and CS,, 1.1405. 
Skewness coefficients, LN2, W2, and GA2 (at variance 0.3957) are: 
2.166, 0.952, and 1.273, respectively. 

Case 1.—Results obtained using CS, for all distributions: this case 
for the four distributions. result, LN, and will enter 
the negative range and will bounded upwards. Since closer 
the case will have better acceptable bounds. The lower 
bounds LN, and are —1.387, —0.081, and —0.556, respectively. 
The upper bound exp (c) 4.306. The effect relatively lower 
upper bound reflected its low values compared the other 
three distributions when 100 (Table 10). The negative areas LN, 
and are 2.62%, 0.86%, and 2.25%, respectively. Although has better 
bounds and lower negative area, its quantiles for ALE are inherently low 
compared and this case, the choice distribution makes 
difference the final results for 100 yr. The quantiles differ from 
2%-8.5% increased from 100 yr-500 yr. However, LP, because 
its low upper bound, and owing its inherently low quantiles are considered 
unsuitable this particular instance, the next best choice but also 
gives comparative results. 

Case 2.—Results obtained using the new unbiased estimates skewness 
for each distribution: this case also, the estimates skewness are less than 
for the four distributions, but are greater than result, the attributes 
the distributions are improved compared Case The lower bounds 
LN, and are —1.030, —0.005, and —0.344, respectively. The upper 
bound exp (c) 7.367. The negative areas LN, and 
are 1.83%, 0.014%, and 1.26%, respectively. Again, will have the best 
attributes, but its quantiles are inherently low (Table 10). The next best choice 
but and also give comparative results (maximum error for 

The foregoing examples show that use the new unbiased estimates 
skewness may not always improve the statistical attributes distribution. 
This may because such distributions are not suitable for analyzing particular 
hydrologic sample. Nevertheless, view the new knowledge regarding the 
behavior sample skewness, use the new unbiased estimates skewness 
preferable place the conventional CS,. 

two instances (Case Sample and Case Sample 2), the 
has seemingly low upper bounds 11.450 and 7.367, respectively), but 
its upper quantiles compare well with the best distribution identified the 
criteria this paper. This suggests area for further investigation, i.e., whether 
can emerge the single most practical distribution for application, which 
may always give results comparable the best alternative distribution (especially 
when the method used for parameter estimation). Note that when 
has upper bound, its lower tail starts the zero value which condition 
most suitable for analyzing the ASE. 

computer program available fit data the four distributions. The 
program evaluates different results presented the foregoing examples including 
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Summary ano 


The lognormal, Weibull, Pearson type and log Pearson type each 
three-parameter distribution, were evaluated generalized fashion terms 
the dimensionless variate which has population mean 
unity. The bounds the distributions, areas the portion distributions 
the negative region variate when they enter such regions, and the differences 
important quantiles among the four distributions are presented the basis 
and y,. Since and y,, this study allows understanding 
the differences existing among distributions and their properties the basis 
statistical parameters (mean, variance, and skewness coefficient) data 
sample. The conclusions from this study are follows: 


The four distributions studied become less applicable for hydrologic fre- 
quency analysis they deviate more and more from their two-parameter 
counterparts (LN2 the case LP). 

given the lower bound negative when less than the 
two-parameter probability distribution value and positive when greater 
than the two-parameter probability distribution value for and when 
positive, and for general. The magnitudes the lower bound 
(positive negative), and the negative area the three distributions increase 
with the deviation from the two-parameter probability distribution value. 

given y,, the lower bound negative when greater than the 
two-parameter probability distribution value and positive when less 
than the two-parameter probability distribution value for and when 
positive, and for general. The magnitudes the lower bound 
(positive negative) and the negative area the three distributions increase 
with the deviation from the two-parameter probability distribution value. 

given the bounded upwards when less than that 
LN2 and the bound value decreases with y,. When greater than 
that LN2, bounded the lower tail and the bound value increases 
with y,. 

given positive y,, bounded upwards when greater than 
that LN2 and the bound value decreases with the increase When 
less than that LN2, bounded the lower tail and the bound 
value increases with the decrease 

When negative, LN, and are bounded upwards and the 
bound values appear generally low. 

will serve rule thumb determine which three-parameter probability 
distribution has better properties for application hydrology; general, the 
three-parameter probability distribution whose two-parameter counterpart (LN2 
case LP) has its close will have better applicable properties. 

For some combinations and the upper quantiles some all 
the four distributions are practically equal for given Thus, the choice 
distribution not critical under those conditions for flood flow frequency 
analysis. 

For some samples, the may have seemingly low upper bound, but 
its upper quantiles may compare well with the best alternative distribution 
identifiable the criteria this paper. 
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The following symbols are used this paper: 


parameters distributions; 


sample skewness coefficient; 
probability; 
F(x) distribution function; 
f(x) probability density function; 
GA2 gamma distribution; 
K,k dimensionless variate; 
constant; 
LN2 two-parameter lognormal distribution; 
log Pearson type distribution; 
sample size; 


Pearson type distribution; 
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w 
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Subscripts 
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sample variance; 

return period; 

three-parameter Weibull distribution; 
two-parameter Weibull distribution; 
random variable; 

sample mean 

gamma function; 

population skewness coefficient; 
population coefficient variation; 
population mean; 

population standard deviation; and 
population variance. 


biased estimate; 

ith value (for variates); 

two-parameter probability distribution; and 
unbiased estimate (conventional). 
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Notes 


provide place within ASCE for publication technical ideas that have not advanced, 
yet, the point where they warrant publication Proceedings paper Journal, the 
publication Technical Notes was authorized the Board Direction October 16-18, 
1967, under the following guidelines: 


original manuscript and two copies are submitted the Manager Technical 
and Professional Publications, ASCE, 345 East 47th Street, New York, N.Y., 10017, along 
with request the author that considered Technical Note. 

The two copies will sent appropriate Technical Division Council for review. 

the Division Council approves the contribution for publication, shall returned 
Society Headquarters with appropriate comments. 

The technical publications staff will prepare the material for use the earliest possible 
issue the Journal, after proper coordination with the author. 

Each Technical Note not exceed pages the Journal. approximation, 
each full manuscript page text, tables, figures the equivalent one-half Journal 
page. 

The Technical Notes will grouped special section each Journal. 

Information retrieval abstracts and key words will unnecessary for Technical Notes. 

The final date which Discussion should reach the Society given footnote 
with each Technical Note. 

Technical Notes will not included Transactions. 

10. Technical Notes will included annual and cumulative subject and author 
indexes. 


The manuscripts for Technical Notes must meet the following requirements: 


Titles must have length not exceeding characters and spaces. 

The author’s full name, Society membership grade, and footnote reference stating present 
employment must appear the first page the manuscript. Authors need not Society 
members. 

The manuscript submitted original copy (with two duplicates) that typed 
double-spaced one side (220-mm) 11-in. (280-mm) white bond paper. 

All mathematics must typewritten and special symbols must properly identified. 
The letter symbols used must defined where they first appear, figures text, and arranged 
alphabetically Appendix.—Notation. 

Standard definitions and symbols must used. Reference must made the lists 
published the American National Standards Institute and the Authors’ Guide the 
Publications ASCE. 

Tables must typed double-spaced (an original ribbon copy and two duplicate copies) 
one side 8-1/2-in. (220-mm) (280-mm) paper. explanation each table 
must appear the text. 

Figures must drawn black ink one side (220-mm) 11-in. (280-mm) 
paper. Because figures will reproduced with width between in. (76 mm) 4-1/2 
in. (110 mm), the lettering must large enough legible this width. Photographs must 
submitted glossy prints. Explanations and descriptions must made within the text 
for each figure. 

References cited text must typed the end the Technical Note alphabetical 
order Appendix.—References. 

Dual units, U.S. Customary followed (International System) units parentheses, 
should used throughout the paper. 
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TECHNICAL NOTES 


COMPARISON Two SURFACE 
HEAT EXCHANGE 


Richard Noble’ 


Mathematical modeling water temperature has received considerable atten- 
tion over the years (1,2,3,5,9,11,12,13). has been observed that not possible 
develop analytic solution for the water temperature many cases due 
the nonlinearity some terms the expression for the net heat flux 
the air-water interface. One example the dependence shown the back 
radiation term. Also, the incoming short-wave radiation term subject rapid 
and unpredictable changes due cloud cover and other atmospheric changes. 
Therefore, methods have been developed linearize the net heat flux expression 
simplify the calculations. 

The concept the equilibrium temperature was developed linearize the 
net heat flux expression. The equilibrium temperature was defined the 
temperature which the water surface would attain when the net heat flux became 
zero. The net heat flux was then expressed the product surface exchange 
coefficient and the difference between the actual water surface temperature 
and the equilibrium temperature. Gameson, Gibbs, and Barrett (7) used this 
concept model the temperature heated river. They used simplified 
analysis with empirical expression relating air temperature and equilibrium 
temperature. Velz and Gannon (14) applied this method heat balance for 
ponds and streams. Edinger, Duttweiler, and Geyer (5) used the concept 
equilibrium temperature detailed linearization the net heat flux term. 
Dingman (4) developed simplified method for computing the equilibrium 
temperature developing curves heat exchange rate (excluding shortwave 
radiation) versus the difference between the water surface temperature and 
the air temperature. The slope and intercept these curves could then 
used predict the equilibrium temperature. 

The equilibrium temperature method has been applied predictive models 
for water temperature. Edinger (6) used this concept model vertical temperature 
structure static systems. Brocard and Harleman (1) developed one-dimensional 
model for water temperatures rivers unsteady flow. 

Prof., Chemical Engrg. Dept., Box 3295, Univ. Station, Univ. Wyoming, 
Laramie, Wyo. 82071. 

Note.—Discussion open until August 1981. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication May 1980. 
This paper part the Journal the Hydraulics Division, Proceedings the American 
Society Civil Engineers, Vol. 107, No. HY3, March, 1981. ISSN 0044- 
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Yotsukura, Jackman, and Faust (15) developed linearization method for 
the net heat flux term expanding each term Taylor series about 
undefined base temperature and neglecting higher than first-order terms. This 
allows one write the net heat flux the sum zero order and first 
order term with respect the water surface temperature. 

Paily, Macagno, and Kennedy (12) used the aforementioned linearization 
method develop analytical solution the dissipation excess heat from 
rivers under subfreezing conditions. The writer (10) developed analytic solution 
river water temperature based the aforementioned linear approximation 
the net heat flux. 

The objective this paper show that relationship can developed 
which relates the parameters the equilibrium temperature model the 
parameters the Taylor Series expansion model. This allows one use either 
method for any model considered interchanging the parameters shown. 
will further shown that the equilibrium temperature model special 
case the Taylor Series expansion model. 

should noted that either linearization scheme mentioned approxima- 
tion and strictly valid only during steady meteorological conditions. 


Description Net Heat 


The net heat exchange air-water interface can expressed 
H=H,-H,- H,- H. H,. 


which sum the short-wave and long-wave atmospheric radiation 


absorbed the water; radiation from the water surface; heat 
transfer due the latent heat vaporization (evaporation); heat transfer 
due conduction the air-water interface; and H,. sensible heat transfer 
due the advection water through the air-water interface. 

The parameter normally very small and neglected practice. The 
expression for then becomes 


Substituting Eqs. and into Eq. one obtains 
Applying the method Yotsukura, Jackman, and Faust (15) Eq. one 


expands Taylor series about undefined temperature and neglects 
higher than first order terms: 
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which net heat exchange term evaluated and 


de, 


Substituting Eqs. and into Eq. 


note that the values and depend the value chosen. 


Edinger, Duttweiler, and Geyer (5) developed linearization Eq. using 
the concept equilibrium temperature. Initially, binomial expansion 
was calculated. All terms higher than second order were neglected: 


Next, the rate change saturation water vapor pressure with respect 
water temperature approximated 


should noted that this not good assumption. Equation can 
rearranged yield 


Substituting Eqs. and into Eq. 


The equilibrium temperature defined the surface water temperature 
which would attained when 


Subtracting Eq. from Eq. 


Neglecting the nonlinear terms 
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Neglecting the second order term Eq. 16, one can solve for 


Comparison Two 
First, note from Eq. that assumed that 
they can used interchangeably. Next, calculate and using 


Substituting the aforementioned into Eq. 


Therefore, the parameters the equilibrium temperature model can evaluated 
terms the Taylor Series expansion model setting 

Therefore, the choice usually not optimum one. Use 
the Taylor Series expansion model with better choice e.g., the natural 
temperature the water unaffected artificial heat, suggested Yotsukura 
al. (15) should yield better results. 

example for comparison, let use the following situation, 279 
0.61 mb/K; and 1.456 mb/K. Using Eq. —549 W/m’. From 
Eqs. and 11, 285 and 41.7 W/m’ the Taylor Series 
expansion method. addition, 40.6 W/m’-K and when 
using Eqs. and for the equilibrium temperature method. expected, 


ve 
K=0, 
Comparing and 24: 
a, 
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the Taylor Series expansion model (Eq. 11) predicts that —549 
20°C. The equilibrium temperature model (Eq. 20) predicts —605 
W/m’, error 9.3%. 

should noted that Eqs. and can equated and general relationship 
developed between the parameters the two models. 


Two methods for linearizing the net heat flux air-water interface were 
considered. The derivation the linearized heat flux was developed for each 
method and compared. was determined that the parameters each method 
could simply related when was set zero the Taylor Series expansion 
method. 

Since setting equal zero normally too low value, recommended 
that the Taylor Series expansion method used with better choice 
such the natural water temperature. 
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RETURN PERIOD FOR MEAN ANNUAL 
EVENT 


INTRODUCTION 


Engineering hydrologists sometimes wish assign return period (7) 
the mean annual hydrologic event. The value the mean annual event usually 
the arithmetic average (X) annual series data, such annual flood flows, 
annual low flows, etc. The value not constant for this event (X), 
but depends upon the probability distribution (PD) assumed for the data sample. 
represents the nonexceedance probability for based the assumed 
PD, for given 1/(1 data are annual largest events 
(ALE), such flood flows and 1/F data are annual smallest events 
(ASE), such low flows. For the normal and Gumbel distributions which 
have constant skewness coefficient (y), also constant; for the normal 
distribution for both ALE and ASE, and for the Gumbel 
2.33 when applied ALE. For other PDs, the value depends 
the statistical parameters the random variable. This paper presents 
values for some commonly used PDs generalized fashion. 


Given data can made dimensionless the transform which 
represents dimensionless variate corresponding the ith data item X,. 
The mean dimensionless data 1.0. Different PDs can evaluated 
terms assuming the population mean the random variable 


Engrg. Div., Water Resources Dept., St. Johns River Water Management 
Dist., P.O. Box 1429, Palatka, Fla. 32077. 

Note.—Discussion open until August 1981. extend the closing date one month, 
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This paper part the Journal the Hydraulics Division, Proceedings the American 
Society Civil Engineers, Vol. 107, No. HY3, March, 1981. ISSN 0044- 


HY3 TECHNICAL NOTES 367 


coefficient variation The coefficient skewness being dimensionless, 
common for and The two- and three-parameter lognormal (LN), gamma 
(GA), and Weibull (WB) distributions and log Pearson Type (LP) distribution 
were evaluated for wide range and The value was determined 
1.0 and values were then calculated for ALE and ASE. Note 


TABLE 1.—Values for Two Parameter Probability Distributions 


Annual Largest Events Annual Smallest Events 


Skewness 
coefficient 


that the three-parameter equivalent the Pearson Type distribution. 
For brevity, equations PDs are presented herein. Different relations for 
the four distributions are available Ref. (For LN, GA, and WB, the equations 
three-parameter distributions can converted into the equations two- 
parameter distributions setting the location parameter, 0.) 


Weibull Lognormal Weibull 
(1) (2) (4) (5) (7) 
0.05 2.19 2.13 1.92 1.84 1.89 2.09 
0.10 2.28 2.18 2.00 1.78 1.84 2.00 
0.20 2.41 2.27 2.12 1.71 1.79 1.89 
0.30 2.51 2.34 2.22 1.66 1.75 1.82 
0.40 2.59 2.40 2.31 1.63 1.71 1.76 
0.50 2.67 2.46 2.39 1.60 1.68 1.72 
0.60 2.73 2.52 2.46 1.58 1.66 1.68 
0.70 2.79 2.57 2.53 1.56 1.64 1.65 
0.80 2.85 2.62 2.60 1.54 1.62 1.63 
0.90 2.90 2.67 2.66 1.53 1.60 1.60 
1.00 2.95 2.72 2.72 1.51 1.58 1.58 
2.00 3.33 3.15 3.20 1.43 1.46 1.46 
3.00 3.60 3.54 3.56 1.39 1.39 1.39 
4.00 3.80 3.90 3.86 1.36 1.34 1.35 
5.00 3.97 4.25 4.11 1.34 1.31 1.32 
TABLE 2.—Values for Three Parameter Probability Distributions 
Annual Largest Events Annual Events 
Lognormal Weibull Lognormal Weibull 
(1) (2) (4) (5) (7) 
0.50 2.14 2.14 2.15 1.88 1.88 1.87 
1.00 2.29 2.31 2.33 1.78 1.77 1.75 
1.50 2.43 2.50 2.52 1.70 1.67 1.66 
2.00 2.55 2.72 2.72 1.64 1.58 1.58 
3.00 2.78 3.25 3.08 1.56 1.44 1.48 
4.00 2.95 3.90 3.41 1.51 1.34 1.41 
5.00 3.10 4.66 3.70 1.48 1.27 1.37 
7.00 3.33 6.44 4.20 1.43 1.18 1.31 
9.00 3.50 8.54 4.62 1.40 1.13 1.28 
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For two-parameter LN, GA, and WB, values are tabulated function 
(Table 1). For the three-parameter versions the preceding three 
PDs, function only (Table 2). For LP, depends both 


TABLE 3.—Values for Log Pearson Type Distribution 


Skewness Coefficient 


Note: U-shaped distribution. 


(or and (Table 3). For ALE, increases with and both, 
but for ASE decreases the latter are increased (with some exceptions 
the case LP). For practical application Tables 1-3 one may first determine 
the sample statistics mean variance (S?), and the skewness coefficient 


(1) (7) (8) 
(a) Annual Largest Events 
0.05 2.14 2.28 2.41 2.54 2.78 3.19 3.78 
0.10 2.15 2.28 2.41 2.52 2.71 3.01 3.39 
0.20 2.17 2.31 2.42 2.52 2.68 2.90 3.14 
0.30 2.19 2.33 2.45 2.54 2.68 2.87 3.06 
0.40 2.21 2.36 2.48 2.57 2.70 2.87 3.03 
0.50 2.23 2.39 2.51 2.60 2.73 2.88 3.03 
0.60 2.26 2.43 2.55 2.63 2.75 2.90 3.03 
0.70 2.46 2.58 2.67 2.79 2.92 3.04 
0.80 2.50 2.62 2.70 2.82 2.95 3.06 
0.90 2.53 2.66 2.74 2.85 2.97 3.08 
1.00 2.57 2.70 2.78 2.89 3.00 3.10 
2.00 3.21 3.26 3.31 3.35 
3.00 3.67 3.62 3.60 
5.00 4.69 4.29 4.09 
Annual Smallest Events 
0.05 1.88 1.78 1.71 1.65 1.56 1.46 1.36 
1.87 1.78 1.71 1.66 1.58 1.50 1.42 
0.20 1.86 1.77 1.70 1.66 1.60 1.53 1.47 
0.30 1.84 1.75 1.69 1.65 1.59 1.54 1.49 
0.40 1.83 1.73 1.68 1.64 1.59 1.54 1.49 
0.50 1.81 1.72 1.66 1.63 1.58 1.53 1.49 
0.60 1.80 1.70 1.65 1.61 1.57 1.53 1.49 
0.70 1.68 1.63 1.60 1.56 1.52 1.49 
0.80 1.67 1.62 1.59 1.55 1.51 1.48 
0.90 1.65 1.60 1.57 1.54 1.51 1.48 
1.00 1.64 1.59 1.56 1.53 1.50 1.48 
2.00 1.45 1.44 1.43 1.43 
3.00 1.37 1.38 1.38 
4.00 1.32 1.34 1.35 
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(CS). Then, the sample coefficient variation given by, S,/X. With 
CV’ and one may enter Tables 1-3 and obtain approximate 
value for When negative, values ALE and ASE Table 
represent the values ASE and ALE, respectively, for and GA. 

the case LP, for some combinations and the parameter values 
not permit evaluation values, and for some combinations the form 
shaped. both instances values are not shown Table 

application the results presented this paper follows: Different 
regulatory agencies apply certain criteria allowing development within the 
riverine floodplains. Assume that one such agencies does not permit any 
development activity within the mean annual floodplain rivers its jurisdiction. 
The structures built before the enactment such regulations escape the law. 
Suppose river within the jurisdiction the agency has 1.0 and 
2.0 for its floods. Property located within the mean annual floodplain will 
inundated, the average, once 2.95 floods are distributed the 
two-parameter LN, but once 2.55 distributed the three-parameter 
LN. insurance agency accepting flood insurance may charge more insure 


such property assumes the three-parameter determining its premium 
rates. 


Rao, V., Probability Distributions,’’ Journal the Hydraulics 
Division, ASCE, Vol. 107, No. HY3, Proc. Paper 16124, Mar., 1981, pp. 339-358. 
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Discussions 


Discussions may submitted any Proceedings paper technical note published any 
Journal any paper presented any Specialty Conference other meeting, the Proceedings 
which have been published ASCE. Discussion paper/technical note open 
anyone who has significant comments questions regarding the content the 
note. Discussions are accepted for period months following the date publication 
paper/technical note and they should sent the Manager Technical and Professional 
Publications, ASCE, 345 East 47th Street, New York, N.Y. 10017. The discussion period may 
extended written request from discusser. 

The original and three copies the Discussion should submitted (220-mm) 
(280-mm) white bond paper, typed double-spaced with wide margins. The length 
Discussion restricted two Journal pages (about four typewritten double-spaced pages 
manuscript including figures and tables); the editors will delete matter extraneous the 
subject under discussion. Discussion over two pages long will returned for shortening. 
All Discussions will reviewed the editors and the Division’s Council’s Publications 
Committees. some cases, Discussions will returned discussers for rewriting, they 
may encouraged submit paper technical note rather than Discussion. 

Standards for Discussions are the same those for Proceedings Papers. Discussion 
subject rejection contains matter readily found elsewhere, advocates special interests, 
carelessly prepared, controverts established fact, purely speculative, introduces personalities, 
foreign the purposes the Society. All Discussions should written the third 
person, and the discusser should use the term when referring himself. The 

Discussions have specific format. The title the original paper/technical note appears 
the top the first page with superscript that corresponds footnote indicating the 
month, year, author(s), and number the original paper/technical note. The discusser’s full 
name should indicated below the title (see Discussions herein example) together with 
his ASCE membership grade (if applicable). 

The discusser’s title, company affiliation, and business address should appear the first 
page the manuscript, along with the Proceedings paper number the original 
note, the date and name the Journal which appeared, and the original author’s name. 

Note that the discusser’s identification footnote should follow consecutively from the original 
paper/technical note. the paper/technical note under discussion contained footnote numbers 
and the first Discussion would begin with footnote and subsequent Discussions would 
continue sequence. 

Figures supplied the discusser should designated letters, starting with This also 
applies separately tables and references. referring figure, table, reference that 
appeared the original paper/technical note use the same number used the original. 

suggested that potential discussers request copy the ASCE Authors’ Guide 
the Publications ASCE for more detailed information preparation and submission 
manuscripts. 
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The writers wish thank Almquist and Holley for their interests the 
paper and for their comments. 

was recognised the writers that the transverse distribution velocity 
always plays important role the differential advection natural rivers. 
However, the rectangular flume which the reported experiments were 
performed, where two-dimensional dead zones were present the form 
strips across the flow, the following points are noted: (1) The vertical velocity 
distribution varied very little across the flow; and (2) tests carried out this 
flume with six probes mounted across the flow equal intervals demonstrated 
little difference the statistical parameters calculated for these different 
locations. 

The writers were therefore satisfied that the flow closely approximated 
two-dimensional situation, and, for this case, that comparison between physical 
measurements made one point flow with cross-sectionally averaged 
computations the numerical model was justified. had also been found 
that there was very little difference the values parameters calculated from 
higher moments (i.e., varience and skew): (1) Where these were computed 
employing only the moments from the flow zone elements comparable 
the probe; and (2) where the figures represented the whole cross 
section. 

The foregoing situation applied only the two-dimensional experiment with 
dead zones. the case where was examined, the velocity 
gradient across the flow was greater shown Fig. 13. Similar transverse 
variation concentrations was then found those reported the discussers 
Fig. 12. The writers’ observation this transverse variation led more 
refined consideration the physical relationship the probe positions and 
dimensions the three-dimensional model elements discussed Ref. 27. 

was never intended that results from the two-dimensional (vertical) model 
should compared directly natural streams they were produced 
demonstrate the importance the dead zone concept and provide physical 
understanding the mechanism related the delay the Taylor period caused 
dead zones. 

Following velocity measurements made water race (27) the writers 
calculated the value Eq. for the side dead zones the average 
the values those flow elements horizontally adjacent and corresponding 
closely possible width the flow depth. Physical observation suggested 
this was appropriate. While the writers agree with the discussers’ observation 

1979, Eric Valentine and Ian Wood (Proc. Paper 14790). 
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that transverse models the dead zone mass exchange controlled more 
parts the velocity distribution, the lower side dead zone exchange velocity 
led longer side dead zone retention times. This emphasizes the importance 
the mechanism the dispersion process streams. 

The short comings the assumption are apparent but 
progressing the three-dimensional case including side dead zones practical 
method comparing experimental with numerical computations was required. 
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FIG. 13.—Smooth Bed: Velocity Gradient Across Flow 


agreed that the ‘‘frozen assumption appears met the 
two-dimensional experiments only for the higher values The spatial domain 
experiments were performed investigation the errors inherent employing 
the cloud’’ conversion where dead zones are present. The errors were 
small and were considered justify its use until the experimental and numerical 
results could compared directly the same domain. 

answer the final points raised the discussers: 


30 30 
lume side Flume 
0 50 
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The difficulties making accurate depth measurements where very 
small prompted the comparison experimental data with numerical runs 
employing mm. Some small standing waves increased this problem 
especially for the smallest depths. The velocity distribution the laboratory 
flume was thoroughly measured for those depths where miniature flow meter 
could accommodated. Values were then computed. 

The increase the values (Table with decreasing y,/d are due 
the use the depth (y,,) the nondimensionalising length scale. This reverses 
the progression values which ranged from 4.65 (y, 
(y, mm). The duration for these extreme cases were approximated 
sec and sec, respectively. This would appear agreement with the 
discussers’ observations transverse model slight decrease the duration 
the initial period for decreasing 


27. Valentine, M., Rough Rectangular Journal Hydraulics 
Division, ASCE, Vol. 105, No. HY12, Proc. Paper 15088, Dec., 1979, pp. 1537-1553. 


Errata.—The following correction should made the original paper: 


Page 1010, Table Col. heading: Should read ‘‘y,, instead 


The writer expresses his appreciation the two discussers for their interest 
the paper and their thoughtful comments and questions. One question raised 
Talukder was; imperative make model study for evey pumping 
the discusser pointed out, even when pumping pit built which 
appears identical one that was modeled and operated successfully, 
one may still encounter operational problems. Also model cannot guarantee 
problem-free prototype. The decision whether not conduct model 
study for each installation complex one. However, the final decision appears 
the writer depend primarily the level peace mind the owners 
wish have. somewhat like buying insurance policy. thorough model 
study minimizes the chances difficulties. Then difficulties are encountered 

1979, Paul Tullis (Proc. Paper 14812). 
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the designers and owners know that they took all normal precautions. 
more direct answer the question, the writer recommends that with present 
technology and the lack correlation and generalization model study data, 
model study performed each major pumping pit installation. This will 
only changed when the results the many model studies are correlated 
and additional research funded provide more fundamental information 
the problem. 

Talukder’s second question related who should initiate the model study. 
The discusser answered his own question quite well. The two parties who are 
expected provide expert advice the owners are the consulting engineer 
and the pump manufacturer. Either are position recommend model 
study. However, the involvement the pump manufacturer may come too 
late. may not enter the picture until after the civil works are completely 
designed. is, therefore, suggested that the major responsibility for considering 
the modeling need rests the consulting /design engineer. 

The third comment Talukder was how reduce the need for many 
models. The writer and discusser again agree general solution. Funds must 
made available allow indepth comparison many recent studies. Basic 
research also needed the physics the problem and look scale 
effects problem. The writer has attempted obtain Federal funding for such 
study without success. Unless the government cooperative effort 
pump companies and engineering firms provide such funding, the need modeling 
each installation will continue. 

The second discusser, Spurr, emphasized the importance vorticity and 
the influence turbulence. Vorticity particularly nasty problem deal 
with the model since there not satisfactory means direct measurement. 
Intrusive devices upset individual vortex and since the vortex not stationary, 
nonintrusive devices cannot used. qualitative estimate from visual observa- 
tions still the accepted means classifying individual vortices. 

Proper modeling the turbulence level definitely important evaluating 
the susceptibility design generate and sustain vortices. The largest single 
problem, viewed the writer, that most model studies are performed 
before the prototype built and operated. Consequently, not possible 
obtain prototype turbulence measurements guide the operation the model. 
There also the problem quantifying the dependence the vorticity 
turbulence and how scale this relationship down the model. would 


helpful have more model-prototype studies evaluate the influence 
turbulence. 


DISCUSSION 


and John Kennedy,’ ASCE 


The writers concur with Glover’s view that physical models have been 
invaluable tool river-engineering practice for many years, and that they will 
continue play central role. Nevertheless, the writers still are the opinion 
that there substitute for good prototype observations and data, and under 
optimal conditions prototype and model investigations will proceed consort. 
sure, many riverine problems cannot solved with field data alone, 
although some can. the same token, model data produced without calibration 
and verification using prototype results must viewed with high degree 
circumspect. 

Glover properly questions the reliability the bed-load-transport data obtained 
with the Helley-Smith sampler. the time the reported field data were obtained, 
relatively little work had been done establish the accuracy results obtained 
with this device. the intervening years the present, however, data have 
been published which largely remedy this deficiency. his detailed field 
calibration the Helley-Smith sampler the East Fork River, Wyo., Emmett 
(17) found that the sampler has near-perfect sediment trapping efficiency 
for the sediment particle sizes between 0.50 and mm. Most the bed-load 
sediment sizes the writers’ study areas fall this range. field study similar 
the writers’ was conducted 1978 the same Pool reaches, calibrate 
several one-dimensional numerical sediment models (18). The bed-load data 
obtained 1978 are plotted Fig. against the mean flow velocity, together 
with the 1976 data included the writers’ paper. Note that the measured bed-load 
discharge given units tons/ft/day. seen the figure, the bed-load- 
discharge data vary over nearly two log cycles. The correlation between two 
variables seems astonishingly good, considering the wide variations 
lateral distributions the bed-load discharge natural rivers, and the two 
sets data are good agreement. While this agreement does not conform 
the accuracy the Helley-Smith data, does lend credibility them. should 
added that elaborate laboratory testing the Helley-Smith samplers has 
also been undertaken the St. Anthony Falls Hydraulic Laboratory, University 
Minnesota. 

Physical model, numerical, and prototype investigations will all continue 

1979, Tatsuaki Nakato, Richard Baker, and John Kennedy (Proc. 
Paper 14978). 

Scientist, Inst. Hydr. Res., and Adjunct Asst. Prof., Energy Engrg., 
Univ. Iowa, City, lowa 52242. 
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FIG. 12.—Relationship between and 


essential river-engineering tools, and many problems will require application 
all three. Nevertheless, the writers believe that well-known limitations 
physical and numerical modeling techniques and the results they produce will 
make prototype study the preferred approach some instances, especially those 
involving complex, three-dimensional, bed-load-transporting flows. 


17. Emmett, W., Field Calibration the Sediment-Trapping Characteristics 
the Helley-Smith Bedload Sampler,’’ Open-File Report 79-411, United States Geological 

Survey, 1979. 

Vadnal, L., “‘Field Study Sediment Transport Characteristics the Mississippi 

River Near Buzzard thesis presented The University 


City, lowa, 1979, partial fulfillment the requirements for the degree Master 
Science. 
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Errata.—The following correction should made the original paper: 


Page 1379, Table heading: Should read Principal Quantities 


Obtained from Field instead Principal Quantities Obtained 
from Field Data 
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DISCUSSION 


Frequency Estimates ALLUVIAL 


The writer fears that the discusser did not read the paper carefully, and 
the discusser’s conclusions are therefore somewhat too harsh. The paper was 
not meant treatise the development alluvial fans. The excellent 
references Bull are not germane the problem the determination 
flood hazards alluvial fans. discharge usually passes through the 
fan apex and the uppermost segments the fan well-defined, perhaps 
entrenched, not felt the writer correct helpful unless 
and perhaps entrenched’’ are more quantitatively 
defined. 

Many (most) fans are not entrenched, and how well channel defined 
alluvial fan normally the eye the beholder. Where fans are entrenched 
the flood frequency distribution the apex can considered apply through 
the entrenched portion, and the method described the writer applied that 
portion below the lower point entrenchment. 

The discusser makes four major points: They shall answered order. 


flood, calculated for the fan apex cannot have recurrence interval 
significantly less than 100 the intersection point. For the 100-yr flood 
the apex be, say, the 50-yr flood lower point the fan requires 
that two floods that magnitude greater must pass that point. Those two 
floods must pass the apex. Therefore, the 100-yr flood the apex fact, 
50-yr flood, quod erat demonstrandum. 

The development alluvial fans random process. Regularity the 
mean value function process does not require deterministic process. 
The discusser referred the work Price (16) for excellent discussion 
the theories the formation alluvial fans. The writer does not see how 
this multitude contrasting and conflicting theories ‘‘absolutely refute any 
assumption pertaining random development model for alluvial The 
writer does concede that the various geomorphologists attempt refute each 
other. Perhaps the discusser does not understand that 
and ‘‘somewhat are almost synonymous. The name random 
variable derives from the fact tnat its value cannot generally predicted exactly 
advance (15). valid question might whether, fact, should 
that channel caused given flood has equal probability cross 
any point This assumption uniform distribution somewhat 
arbitrary, and quantitative evidence the contrary would appreciated. 

The assumption log-Pearson III distribution the apex operational 


“November, 1979, David Dawdy (Proc. Paper 15001). 
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assumption, stated the writer. That distribution, however, used only 
for convenience interpolation. The choice distribution nonparametric 
distribution could used) not inconsistent with either the statements 
quoted the discusser. 

Inthe conclusions section, the writer states the three assumptions underlying 
the method. Each those assumptions may counterargued against. However, 
rather unfair the discusser quote out context accuse 
the writer aside’’ the problems which concedes and with which 
ends his paper. Suggestions for improvement would appreciated, particularly 
quantitative, and based field evidence germane the problem 


determination flood frequency distributions arbitrary point alluvial 
fan. 


Cramer, H., and Leadbetter, R., Stationary and Related Stochastic Processes, 
John Wiley and Sons, Inc., New York, N.Y., 1967, page 

16. Price, E., Jr., Random-Walk Simulation Model Alluvial Fan 
Technical Report No. University Arizona, 1972. 


Errata 


The following corrections should made the original paper: 


Page 1395, Eq. Should read constant (C,) 


b 


instead constant (C,) 


Page 1395, Eq. Should constant (C,) 


instead constant (C,) 


Page 1400, Fig. Should Eq. and Eq. instead Eq. and Eq. 
respectively. 


“November, 1979, Steven Wright and Roger Wallace (Proc. Paper 15003). 


DISCUSSION 


The authors have presented application the Galerkin finite element 
method important problem applied hydrology, namely, the simulation 
time-dependent rainfall-runoff. However, the authors’ model presently 
formulated includes deficiency serious enough render the model unacceptable 
for practical applications. 

The authors have solved the time-dependent, two-dimensional equations 
momentum and continuity for rainfall-runoff. For the hypothetical problem shown 
Fig. rainfall-runoff from impervious plane, these equations appear 
properly stated and solved. However, for practical problems such the 
Ina Basin example given the paper, the equation continuity stated 
incorrectly. 

The equation continuity, Eq. accounts for three terms: (1) Rainfall 
the watershed; (2) discharge from the watershed; and (3) the rate change 
flow depth the watershed. What the authors have failed take into 
account here are the significant rainfall losses that occur all watersheds. 
These losses can include infiltration, surface depression storage, and interception 
vegetal cover. Quantifications these losses are routinely included the 
most widely-used rainfall-runoff models, such the United States Environmental 
Protection Agency’s Storm Water Management Model (16). For the authors’ 
model properly represent rainfall-runoff, the equation continuity must 
modified include loss term. 

The example problems which the authors presented avoided the issue rainfall 
losses. Their model reproduced the results Schreiber and Bender (11) well, 
but these data were obtained laboratory experiment using small impervious 
plane model watershed. The authors also applied their model the Ina 
Basin the Tenryu River, but rainfall-runoff data from this site were available 
for purposes comparison. 

summary, the authors have developed finite element model rainfall-runoff 
that flawed because does not account for significant rainfall losses such 
infiltration and depression storage. However, what more important here 
that the authors’ results are example disturbing trend water resources 
modeling which increasingly sophisticated mathematical solutions become 
end themselves. Although advances applied mathematics have enabled 
solve some problems hydraulics and hydrology that could not solve 
before, the importance the mathematics must kept proper perspective. 
not true that the more mathematically sophisticated model is, the better 
model is. this paper, the authors have presented impressive example 

1980, Mutsuto Kawahara and Teruyuki Yokoyama (Proc. Paper 15341). 


Camp Dresser and McKee Inc., 7620 Little River Turnpike, Annandale, Va. 
22003. 
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the numerical solution differential equations, but have failed properly 


state the most basic principle all hydraulic and hydrologic models, the principle 
continuity. 


16. Huber, C., al., Water Management Model: User’s Manual, Version 


Report No. U.S. Environmental Protection Agency, Cincinnati, 
Ohio, Mar., 


Discussion Rema and Arun 


The authors are commended for their effort predict analytically the 
turbulent fluctuations uniform flow channel. The writers, however, feel 
that this study mainly reiterates what already known about turbulence 
characteristics and the formulation the prediction model. 

major part the paper devoted towards deriving working model for 
predicting turbulent intensities velocity for the streamwise and crosswise 
directions open channel uniform flow. The basic assumptions, the subsequent 
derivation, the model and the concluding observations about the values the 
model parameters and are exactly similar the one already published 
the first author along with Shen (3). 

The suggested model relates the turbulent intensity u’, with Taking 
the product and equal —0.4 (Table for smooth boundaries, the 
relationship given Eq. can reduced 


1 1/2 
u 


other words, the maximum value streamwise turbulent velocity 
1.58 whereas, Laufer’s (2) experiments, the maximum value ranged 
from 1.8 Similarly, Blinco and Partheniades (11) observed the range 
maximum value 2.6-3.0 This shows that the model perhaps, 
underestimates the maximum turbulent intensity u’. Moreover, this model which 
expressed decreasing function the relative distance cannot estimate 
near the boundary where increases with y/h. The variation with 
depth warrants two different relationships suggested Nelluri and Novak 
(12). 

1980, Ruh-Ming Li, James Schall, and Daryl Simons (Proc. Paper 
15347). 


Asst. Prof. Civ. Engrg., Delhi Coll. Engrg., Kashmere Gate, Delhi-6, India. 
Prof. Civ. Engrg., Coll. Engrg., Kashmere Gate, Delhi-6, India. 
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The authors have estimated taking the hydraulic radius, equal the 
depth flow, This equivalent assuming that the laboratory flumes 
whose data are studied act hydraulically infinitely wide. For example, the 
measurements taken the authors for discharge 0.0056 0.003 
bed slope with approximate average velocity cm/s (Table 
20-cm wide, 1:1 trapezoidal flume results significantly different the 
bed. Approximating equal gives 3.66 cm/s whereas, using 
directly, 3.16 cm/s. Since the turbulent intensities and are directly 
related u,, important know the reasons for using depth flow 
the hydraulic radius even when the channel not wide. 

The author’s observation that y/h 0.15, u'/u function y/h 
only being negligible) variance with the observation Blinco 
and Partheniades (11) who concluded that apparently universal function 
y/v rather than y/h. 

From Table obvious that the product and approximately 
constant but does not imply that either individually constant throughout 
the depth the channel. Taking the value 0.55 tends overestimate 
transverse fluctuating velocity near the bed and underestimate near the 
free surface (Table 2). Moreover, Smutek (13) has shown that not only 
function y/h but also the Reynolds number. Both these factors could 
have been accounted for simple functional relationship between and 
had been derived. 


Blinco, H., and Partheniades, E., Characteristics Free Surface Flows 
over Smooth and Rough Association Hydraulic Research, 
Vol. No. 1971, pp. 43-72. 


Nalluri, C., and Novak, P., Characteristics Smooth Open Channel 


Circular Cross International Association Hydraulic Research, Vol. 
11, No. 1973, pp. 343-368. 


Hydraulics Division, ASCE, Vol. 95, No. HY1, Proc. Paper 6323, Jan., 1969, pp. 
519-523. 


Errata 


The following corrections should made the original paper: 


Page 584, Fig. the ordinate: should read instead 

Page 599, Appendix II, line 10: Should read distance from sill 
instead distance from the sill 


1980, Rangaswami Narayanan and Loizos Schizas (Proc. Paper 15368). 
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Discussion Don Wood,’ ASCE 


The authors present method solve the pipe network head equations 
simultaneously. They have failed cite several important references which 
also develop procedures for solving the head equations simultaneously. Martin 
and Peters (15) 1963 presented technique which utilized Newton’s method 
gradient approximation linearize the nonlinear head equations. 1966 
Marlow, al. (14) formulated the head equations using approximation which 
appears very similar the method described this paper. Shamir and 
Howard (16), Lui (13) and Lemieux (12) all presented similar methods based 
the simultaneous solution the head equations. appears that, considering 
the similarities the approach pipe network analysis, these references should 
have been cited. Certainly this approach has been widely advocated, and only 
the method utilized linearize the nonlinear head equations seems different. 
also appears that the simplified method proposed the authors inferior 
the use gradient methods. This certainly true the method developed 
this writer and cited the paper where the flow equations were linearized 
similar manner [Wood, Charles (10)]. subsequently found that lineariza- 
tion the flow equations using gradient methods has convergence characteristics 
superior those obtained the method proposed Ref. 10. likely that 
the same true the head equations and the method proposed the authors 
has inferior convergence properties. 

The authors also state that convergence problems occur with this method, 
and that method can guaranteed. While true that method 
guaranteed, several methods based the flow equations have excellent conver- 
gence characteristics which consistently produce reliable results. Studies the 
University Kentucky consisted analyzing over 100 pipe networks ranging 
size over 500 pipes. The correct solution was obtained for every case 
average under seven trials using method based the simultaneous 
solution the flow equations. Using method based the simultaneous solution 
the head equations, convergence difficulties were encountered frequently, 
and failure rate nearly 20% occured. pointed out the authors, lines 
carrying low flows caused particular difficulty. This can generalized state 
that lines with small head changes may cause convergence difficulties for this 
method. Unfortunately, this does not mean that these lines can ignored 
the analysis. Often lines connecting storage tanks have very small head changes 
but carry significant flow and have great effect the system. accurate 
steady state analysis the pipe system requires accurate flow calculation 
for these pipes. appears that algorithms using the head equations are unable 

“July, 1980, Lewis Issacs and Kevin Mills (Proc. Paper 15559). 


Prof. Civ. Engrg., Coll. Engrg., Dept. Civ. Engrg., Univ. Kentucky, Lexington, 
Ky. 40506. 
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this consistently. Others have reported convergence problems associated 
with these methods. Shamir and Howard (16) noted oscillations and reported 
that solution may not obtained. Lui (13) noted that divergence was possible 
and Collins and Kennington (11) presented some data which documented 
convergence problems. This writer has concluded that algorithms are available 
for pipe network analysis which are reliable, but these are based solving 
the flow equations (or flow correction equations). The method presented 
the authors, which based head equations, will unreliable certain 
cases. 

The authors’ method for computing flow rates innovation which 
apparently improves convergence. This may significant improvement 
can shown generally and eliminate significant number 
convergence problems. The authors also noted that accelerator factors were 
applied and different numerical values were used. was not clear how the 
accelerator was applied, but has been the experience this writer that 
accelerators are limited value, because factor which improves convergence 
one situation may adversely affect another. great deal testing 
would required select appropriate accelerator. 


Collins, A., and Kennington, L., discussion Period Simulation 
Journal the Hydraulics Division, ASCE, Vol. 103, No. HY12, Proc. Paper 13378, 
Dec., 1977, pp. 1496-1500. 

Lemieux, F., Algorithm for Distribution Journal the 
Hydraulics Division, ASCE, Vol. 98, No. HY11, Proc. Paper 9336, Nov., 1972, pp. 
1911-1920. 

Liu, H., Numerical Analysis Water Supply Networks Digital 
Proceedings, 13th Congress the International Association for Hydraulic 
Research, Vol. Subject Sept., 1969, pp. 35-42. 

Marlow, A., Hardison, L., Jacobson, H., and Biggs, E., Design 
Fluid Networks with Journal the Hydraulics Division, ASCE, 
Vol. 92, No. HY4, Proc. Paper 4866, July, 1966, pp. 43-61. 

Martin, W., and Peters, G., Application Newton’s Method Network 
Analysis Digital Journal the Institute Water Engineers, Vol. 
17, 1963, pp. 115-129. 

Shamir, U., and Howard, D., Distribution Systems Journal 
the Hydraulics Division, ASCE, Vol. 94, No. HY1, Proc. Paper 5758, Jan., 1968, 
pp. 219-234. 
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VISUALIZATION SEPARATION OVER SAND 
Discussion Peter and Lam 


The authors are commended for demonstrating the use simple 
and relatively inexpensive technique measure flow separation lengths over 
dunes. Such method step forward providing information for the 
development semi-analytical models for determining dune length such 
that proposed (14). However, the writers feel that the analysis 
the data not entirely correct. addition, the results obtained the authors 
differ from those obtained from recent experiments (paper preparation) 
the writers. 

general for fixed bed, the flow separation length must expressed 
function the following independent variables (using the author’s notation) 


the dunes; length the dunes; median grain size roughness 
the dunes (for uniform grain size distribution); average velocity the 
flow; average depth the flow; density the fluid; viscosity 


the fluid; and acceleration due gravity. Eq. can written 
dimensionless form 


Yalin (13) has shown that for the boundary Reynolds number 
shear velocity), the sand waves are dunes and under these conditions 
the dune length given A/H 27. Accepting the fact that A/H will 
more less constant, the term can eliminated from Eq. 11. general 
and both have effect However, this case, the authors 
have held A/A constant 0.06. The writers have found that this range 


steepness has effect but significant for smaller values A/A. 
Therefore, for the authors’ data 


Froud number given U/V gH. 


1980, Emin Karahan and Allan Peterson (Proc. Paper 15652). 

Engr., Environmental Hydr. Section, Hydr. Div., National Water Research 
Inst., Canada Centre for Inland Waters, P.O. Box 5050, Burlington, Ontario, L7R 4A6 
Canada. 

“Head, Environmental Hydr. Section, Hydr. Div., National Water Research Inst., Canada 
Centre for Inland Waters, P.O. Box 5050, Burlington, Ontario, L7R 4A6 Canada. 
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For small values Froude number, properties the dunes and flow over 
the dunes should independent the Froude number. The authors have 
plotted L/A versus (without keeping constant) and suggested that L/A 
increases with decreasing when 0.2. They also plotted the same data 
L/A versus showing increasing with decreasing The writers wish 


100 


FIG. Dimensionless Separation Zone Length, with Average 
Flow Froude Number 


results 


Raudkivi (2) 
authors’ results 


FIG. 7.—Variation Dimensionless Separation Zone Length, L/A, with Dune 
Steepness A/A 


emphasize that what shown Fig. actually the effect and not 
any Froude number effect. This can substantiated the writers recent 
data Fig. For constant A/A and and with 30,000, constant 
was obtained for values between 0.05 and 0.45. Therefore the writers 
are the opinion that L/A independent and that the authors had 


1.0 
A A 
4 
VgH 


388 


MARCH 1981 HY3 


used parameter Fig. they might have arrived the same conclusion. 

The effect A/A can seen Fig. which the writer’s data are 
plotted for values A/A 0.02, 0.033, 0.05 and 0.07, keeping constant 
and having again 30,000. The results show clearly defined relationship 
between L/A and A/A. This shows that, for given function 
A/A only. 

Finally, the writers have found that the value L/A obtained from their 
own measurements lower than that obtained the authors. The writers 
obtained their measurements using dye injection technique which they consider 
more direct method and was found very satisfactory over the 
flow conditions tested. The authors have found that their value L/A agreed 
well with results obtained others using rectangular reversing steps well 
those obtained Raudkivi (15), using ripple bed form. The writers feel 
that comparison between dunes and rectangular reversing steps not justified 
because the totally different geometry the two boundaries. the other 
hand, the ripple bed form used Raudkivi (15) had value A/A 0.075 
and, can seen from Fig. when A/A 0.075, L/A 0.32. Therefore, 
Raudkivi’s data agree with those obtained the writers. The writers therefore 
suggest that the authors’ value L/A 0.32 for their case A/A 0.06 
too high (see Fig. and that this may due very shallow depths that 
the authors may have had use because their small flume. Unfortunately, 
the authors did not include the basic hydraulic data their paper and thus 
this depth effect could not determined. 


14. J., Stationary Progress Report No. 49, Technical 
University Denmark, August, 1979, pp. 3-10. 


15. Raudkivi, Y., Sediment Ripple Journal the Hydraulics 
Division, ASCE, Vol. 89, No. HY6, Proc. Paper 3692, Nov., 1963, pp. 15-33. 


Errata 


The following corrections should made the original paper: 


Page 1848, Eq. 34: Should read 


instead 


a+i/2 +h iJ ij+t 


“November, 1980, Christopher Koutitas and Brian O’Connor (Proc. Paper 15836). 
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